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Abstract 



We consider the hard leptoproduction of a photon off nuclei up to spin-1. As a new result we 
present here the general azimuthal angular dependence of the differential cross section for a spin-1 
target. Its twist-two Fourier coefficients of the interference and squared deeply virtual Compton 
scattering amplitude are evaluated in leading order approximation of perturbation theory in terms 
of generalized parton distributions, while the pure Bethe-Heitler cross section is exactly calculated 
in terms of electromagnetic form factors. Relying on a simple model for the nucleon generalized 
parton distribution H, which describes the existing DVCS data for a proton target, we estimate the 
size of unpolarized cross sections, beam and longitudinal target spin as well as unpolarized charge 
asymmetries for present fixed target experiments with nuclei. These estimates are confronted with 
preliminary HERMES data for deuterium and neon. 
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1 Introduction 

Exclusive two-photon processes in the light-cone dominated region, i.e., in the generalized Bjorken 
limit, are most suitable for the exploration of the partonic content of hadrons. This comes from 
the fact that the dominant contribution to the amplitude of such processes arises from Feynman 
diagrams in which both photons directly couple to one quark line P . The processes of interest are 
different photon-to-meson transition form factors, i.e., 7*7^*^ M, the production of hadron pairs 
by photon fusion and the crossed processes like 'y*H B^y or 7*-*'' if Bl^l~ . The latter class 
might be denoted as deeply virtual Compton scattering (DVGS). It contains different processes: 
DVGS on a hadron target without and with excitation of the final state due to the leptoproduction 
of a real photon [21 El IH El IHl El , the photoproduction [S] or leptoproduction EU] of a lepton 
pair. 

The factorization of short- and long-range dynamics is formally given by the operator product 
expansion (OPE) of the time ordered product of two electromagnetic currents, which has been 
worked out at leading twist-two in next-to-leading order (NLO) and at twist-three level in leading 
order (LO) of perturbation theory (for references see [7]). However, one should be aware that the 
partonic hard-scattering part, i.e., the Wilson coefficients, contains collinear singularities, which 
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are absorbed in the non-perturbative distributions by a factorization procedure, which has been 
proven at twist-two level ^] . 

The non-perturbative distributions are defined in terms of light-ray operators with definite 
twist sandwiched between the corresponding hadronic states. These distributions, depending 
on the two-photon process under consideration, are sensitive to different aspects of hadronic 
physics. Especially, in DVCS one can access the so-called generalized parton distributions (GPDs). 
The second moment of the fiavor singlet GPD is related to the expectation value of the energy 
momentum tensor. Thus, it contains information on the angular orbital momentum fraction of 
the nucleon spin carried by quarks ^Hj. Moreover, in contrast to ordinary parton distributions, 
measurable in inclusive reactions, GPDs carry also information about the parton distribution in 
transverse direction and might so provide us with a hologramatic picture of the nucleon fT^ ITKj 

uniini. 

Although GPDs are accessible in the leptoproduction of mesons (TH], DVCS is the theoretically 
cleanest tool to probe the partonic content of the nucleon on the level of amplitudes. The DVCS 
process in the leptoproduction of a photon on a proton target has been measured by the HI col- 
laboration jinj in the small xb region (see also |20j), due to the single beam spin asymmetry by the 
HERMES j2Tj and CLAS [221 collaborations as well as due to the charge asymmetry at HERMES 
j2Hl- Unfortunately, in this DVCS process a deconvolution of GPDs from a measured amplitude is 
practically impossible. Thus, a model for GPDs is required that satisfies the constraints coming 
from first principles [H 121 E] , namely, support properties, sum rules and the reduction to the par- 
ton densities in the forward kinematics. Moreover, different authors derived positivity constraints 
(see Ref. [21] and references therein). Although these derivations are intuitively understandable, 
a number of questions can be raised that are hardly to answer. At LO in perturbation theory 
all experimental data are consistent with an oversimplified GPD model 7 that fulfills the first 
principle and also the positivity constraints, derived in Ref. |i?5| . 

Recently, the measurements of beam spin asymmetries in the leptoproduction of a photon on 
neon and deuteron targets have been reported by the HERMES collaboration |2Z] . Certainly, it is 
appealing to employ DVCS for the investigation of the internal structure of nuclei. Although the 
binding energy is negligibly small compared to the virtuality of the exchanged photon, one might 
expect that DVCS observables are affected by the binding forces. While for spin-0 and -1/2 targets 
one might assume that the nucleus GPDs can be expressed in terms of slightly modified proton 
GPDs, in the case of a spin-1 target new GPDs enter the DVCS observables, in which bound 
state effects are dominantly manifested. An appropriate candidate for studying such effects is the 
deuteron, which has been widely used as a target in lepton-scattering experiments. This nucleus 
has been extensively studied in both deep-inelastic [2H] and elastic [23 IHOl IHIj scattering. From 
the theoretical point of view, it would be desirable to have complementary information, which 
would give us a deeper understanding of the binding forces and could hopefully shed some light 
of its effective description and the fundamental degree of freedom in QCD. 

On the theoretical side it is time to study DVCS on nuclei in more detail. For deuterium 
a parameterization of leading twist-two operator matrix elements in terms of GPDs has been 
proposed in Ref. [S2]- First estimates for the beam spin asymmetry have been presented in Ref. 

based on a drastic simplification of the theoretical prediction, and in Ref. [311 ES] by means 
of a convolution model. In Ref. jSHl it has been argued that nuclei GPDs can deliver information 
about the spatial distribution of the strong forces in terms of the fundamental degrees of freedom 
in QCD. In this paper we complete the twist-two sector for the theoretical predictions for a spin-1 
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target, generalize our oversimplified model, and estimate different asymmetries for fixed target 
experiments. 

The outline of the paper is the following: In Section|21the OPE approach is applied to DVCS on 
a target with arbitrary spin. Employing these general results, in Section |21 we evaluate the Fourier 
coefficients appearing in the cross section for a spin-1 target at leading twist and at leading order 
in perturbation theory. In Section |3] we present oversimplified GPD models for nuclei, expressed 
by the proton ones, where bound state effects for longitudinal degrees of freedom are neglected. 
Relying on qualitative properties of proton GPDs, consistent with the DVCS data for the proton 
target, we estimate in Section El by a rough kinematical approximation beam spin, charge, and 
target spin asymmetries for fixed target experiments at HERMES and JLAB kinematics. We then 
provide numerical results for these asymmetries and the unpolarized cross section. Finally, we 
summarize and give conclusions. Appendices are devoted to the parameterization of the deuteron 
electromagnetic form factors, the results for the Fourier coefficients of a spin-1 target, and the 
target mass corrections. 



2 Azimuthal angular dependence of the cross section 

In this Section we evaluate the differential cross section for the leptoproduction of a photon 

l^{k)A{P{) -^l^{k')k{P,)^{q,) (1) 

off a nucleus target A with atomic mass number A and mass Ma. The goal is to express the 
azimuthal angular harmonics in terms of the electromagnetic current and Compton amplitudes to 
twist-two accuracy at LO in The considerations are valid for any target with arbitrary spin 
content. 

The five-fold differential cross section 

2 



da Oi xp^y 



dxAdyd\A^\d(f)dip 1671^ Q^VTT^ 



T 

i3 



Ma . X 

2xA-^ (2) 



depends on the scaling variable xa = Q^/2Pi-gi, where Q"^ = —qf with qi = k — k', the momentum 
transfer square = (P2 ~ -Pi)^? the photon energy fraction y = Pi ■ qi/Pi ■ k and, in general, 
two azimuthal angles. In the following we refer to the rest frame, shown in Fig. The scaling 
variable xa is related to the Bjorken variable xb by 

xb = ~ Axa , (3) 

2MNEy ' ^ ' 

where Mn is the nucleon mass and E is the lepton beam energy. The photon energy fraction y 
and the Bjorken variable satisfy the well know relation 

where sjv is the center-of-mass energy squared for the lepton scattering off a nucleon. Below it is 
sometimes more convenient to use the scaling variable 

^ Xa Xb 
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Figure 1: The kinematics of the leptoproduction in the target rest frame. The 2;-direction is 
chosen counter-along the three-momentum of the incoming virtual photon. The lepton three- 
momenta form the lepton scattering plane, while the recoiled nucleus and outgoing real photon 
define the nucleus scattering plane. In this reference system the azimuthal angle of the scattered 
lepton is 0; = 0, while the azimuthal angle between the lepton plane and the recoiled nucleus 
momentum is (p^ = (f). When the nucleus is transversely polarized (in this reference frame) 
S± = (0, cos$, sm$, 0), the angle between the spin vector and the scattered nucleus is denoted as 

instead of xa, where A^/ corrections have been neglected. Here t] is called the skewedness 
parameter, giving the longitudinal momentum fraction in the t-channel. As we see, in the DVCS 
kinematics it is simply related to the Bjorken like scaling variable ^. To define the azimuthal angles 
we point the virtual photon three- momentum towards the negative ^-direction, (p = (pi — (p^ is 
the angle between the lepton and nucleus scattering planes and ip = (P — (pN is the difference of 
the azimuthal angle <P of the spin vector 

= (0, cosfPsin©, sinfpsin©, cos©), (6) 

giving the magnetic quantization direction for the initial nucleus, and the azimuthal angle ^at of 
the recoiled nucleus as depicted in Fig. [T] 

We consider this process in the (generalized) Bjorken limit, i.e., ~ Pi ■ qi should be large 
compared to Mj^ and A^. Obviously, increasing the atomic mass number A one will violate the 
condition S> Mj^, since is in reality restricted by the experimental settings. It is instructive 
to consider the situation in deeply inelastic scattering. Here the hadronic tensor is given by the 
absorptive part of the forward scattering amplitude. In the kinematical forward case the target 
mass corrections are given by e^, which is independent on the mass number A. In our kinematics 
the situation is more complex P7j. Fortunately, it is shown in Appendix O that within our GPD 
model, specified below, in the DVCS kinematics the target mass corrections^ possess the same 

^We distinguish between power suppressed corrections that arise from multi-parton correlation GPDs and target 
mass corrections. The former ones contain new dynamical information and not much is known about them. The 
latter ones arise from trace subtraction of twist-two operators. They have been elaborated in Ref. (S7| and are 
given as convolution with so-called double distributions or alternatively with GPDs. 
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A scaling as the twist-two contributions. Thus, they have nearly the same size as for a nucleon 
target and we may employ the OPE approach also for realistic and larger A. 

In the one virtual photon exchange approximation the amplitude T is the sum of the DVCS 
Tdvcs and Bethe-Heitler (BH) Tbh ones, displayed in Fig. |21 

= + ITuvcsl' + 1} , (7) 

with the interference term 

I = 'Tdvcs'Tbh + '^vcs'^H, (8) 

where the recoiled lepton (A'), nucleus (S") as well as photon (A') polarization will usually not be 
observed. Each of these three terms in Eq. (|7j) is given by the contraction of DVCS tensor 

T^,(P, A,g|5,5') = dxe'^-'^{P2,S'\Tj,{x/2)U-x/2)\P,,S) (9) 

or /and electromagnetic current 

J„(P,A|^,^0 = ^(A,'5'|j„(O)|Pi,5), ^ e,7/>7,^/., (10) 

i=u,d,s 

with a corresponding leptonic tensor 

I^VCSp = — (-5'"^)-^DVCS X^^oM (^/3i^)^ (11) 
<ll s' 

I'^Buf = -^^BuZ^Jf^Jl^ (12) 



i=^,L-^^j:u^(T..y^i.c. rzi. (13) 



for e 



for 



Here P = Pi + P2 and q = {qi + q'2)/2 and we summed over the photon (A') and lepton (A') 
polarizations, where the second summation is included in the leptonic tensors. For any nucleus 
(or hadron) target a general parameterization of the DVCS tensor at twist-three level and in LO 



5 



approximation of perturbation theory has been given in Ref. j|38j. To NLO accuracy a new Lorentz 
structure occurs due to the flip of the photon hehcity by two units, which is caused by the so-called 
twist-two gluonic transversity [321 EDI E], describing an angular orbital momentum flip by two 
units, as well as twist-four effects jl2]- In the following we only consider its twist-two part, where 
it is expected that for the typical kinematics of present fixed target experiments the twist-three 
contributions are small for different observables, but not for all of them. In this approximation 
the DVCS tensor is parameterized by a vector V/^ and an axial- vector A'^ part 



Q ■ 

T = -V '^a V - - - V "ie 



-pr 

arqp ' y 



P-q'' 



(14) 



where the projector V^u = g^i, — qi^iq2u/(li ■ q2 ensures current conservation. The leptonic tensors 
are evaluated by means of the standard projector technique and read after summation over the 
helicity of the final state lepton (see Fig. ^ 



-^DVCS 



-^BH 



2 (PA;"^ + k'^k'^' -k-k'gf"" - iX e'^^"^ k^k'^) , 
(-^?,^)Tr r"^(A;, k\ A)^^^^V'^{k, k\ A)^, 

(-^7°^)Trr^^(A:,fc',A)|^i^-^7^r, 



(15) 
(16) 

(17) 



where 



r.,(fc,fc-,A) = "°i^-^^ + "-^^^+^^"° 



(A;- A)2 



{k' + A)2 



;i8) 



is the leptonic part of the BH amplitude. 

Employing the parameterizations (jHT|) . ()33|) and (jHH), the contractions of leptonic and DVCS 
tensors result in the kinematically exact expression for the squared BH term (of course, in tree 
approximation) 



IT. 



-8e« 



BH| 



Pl(0)P2(0)Q^A2 



^\q- J k- + k- J q- f -q- J q- -2k- J k- 

CI 



A2 + 2A;- A)2 + 4(A;- A)^ 
4A2 



J - 



(19) 



iX[{Q^ + Ak- A- A^) egAjjt + A^efc(2g+A)jjt] 
2A2 



while the interference term^ (here for negatively charged lepton scattering) 



X 



2 - 2?/ + ?/2 4e6^ 

(0)^2(0) A2Q4 

A(2 - y)y Ae^^ 
^y^T^i (0)^2(0) A2Q4 



A;" 



y 



q-A[ 



(20) 



+ h.c. 



^Note that the sign in front of the Levi-Civita tensors in Eqs. (20) and (23) of Ref. |S] is erroneously. 
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and the squared DVCS amplitude 

ITDVCsr = 8e^ ^ ~ ^ ^ [q.V, q.Vl + q-A, q.A\) 



y s' 

are expanded with respect to l/Q. In contrast to the squared DVCS amphtude the interference as 
well as the squared BH terms have an additional 0-dependence due to the (scaled) BH propagators 

V, ^ ^^^1^ = -^T^ {J + 2Kcosm , (22) 
y[l + e^) 

(k — 1 



where 



^^(i-.-^)(i^|)-a-)(2-4: 



and 



(23) 

with the plus sign taken for the square root in Eq. ((221) • vanishes at the kinematical boundary 
= A^ijj, determined by the minimal value 



2 _ ^^2{1 - xa) {I - VTT^ + _ Mjxl 

4xA(l-XA)+e2 1 -xa + xaM1/Q2 ' 



-^Ln = Q- JT^T^ zr^T-^ ~q ,t ,^2/^2 > (24) 

as well as at 



In the frame chosen the contractions of leptonic, depending on fc, gi, and A, and DVCS, a 
function of g, P, A and the spin vector S', tensors in Eqs. (fTHHTTjl . yield finite sums of Fourier 
harmonics [HIZI: 

1^™!^ - J^^vS^M ( ^ t [=™ ("^) + »™ - m ] . (25) 



n=l 



Tovcsr = ^ jcr^^ + E [cr^^ cos(n0) + .r^^ sin(n0)] I , (26) 



y'Q 



n=l 

3 



where the + (— ) sign in the interference term stands for the negatively (positively) charged lepton 
beam. In the twist-two sector the (j) independent contribution of the squared DVCS amplitude is 



c, 



DVCS 




{2-2y + y 



-A(2-y)?/ 



q-Vi q-A\ + q-Ai q-V^ 



while the Fourier coefficients of the interference term are 
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'^"^ ^ ^>m I 



X 



+A(2-i/)y 



T ^9A 2k^gAjg-Vi"^ 



(28) 



(29) 



expanded in leading order of Note that the cos(30) and cos(20) harmonics in the interference 
and squared DVCS term, respectively, also appear at twist-two level. Since they arise from the 
gluon transversity, they are suppressed by a^/vr. 



3 Fourier coefficients for a spin-one target 

So far the formalism, presented in the previous section, is rather general and at LO of perturbation 
theory it can be easily extended to the twist-three sector. Introducing an appropriate form factor 
decomposition of the electromagnetic current and the DVCS amplitudes, it can be employed for 
any target. The resulting Fourier coefficients for the spin-0 and -1/2 targets are presented at twist- 
three level in Refs. jHl Ej for a pion^ and proton target, respectively. The results for a nucleus 
target follow from the replacements of kinematical variables, form factors Fj(A^) and (set) of 
GPDs: 



XB^XA, e^e, M-^Ma, F,{A^) ^ Ff^iA^) , (30) 

F = {h, E, H, E, } (x, V, A') ^F^ = {h^, E^, H^, E^, } (x, r/, A^) . 

Let us now derive the Fourier coefficients for a spin-1 target in terms of GPDs. The electro- 
magnetic current 

J, = -el-e,P,Gi + {el-Pei^ + ei-Pe^J G2 - e^- ^^i-P^ G3 (31) 

is given by three form factors G'j(A^) with i = {1, 2, 3}, where ei^ (€2^) denote the three polariza- 
tion vectors for the initial (final) nucleus. The form factors Gj(A^) can be measured due to the 

■^Here the variable xb has be defined in the rest frame with respect to the target mass rather than to the nucleon 
mass. 
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elastic scattering of a lepton on a nucleus. For the deuteron their parameterizations are available 
in the literature, see Ref. 1^01 1^ and references therein. 

At twist-two level the amplitudes Vi and Ai, appearing in the parameterization (fT^ of the 
DVCS tensor, can be decomposed in a complete basis of nine Compton form factors (CFFs) 

^(e, A^, = ^n^n^, Hs, n,, n,, n^, n^, n^, n,] (e, a^, q^). (32) 

For convenience we employed here the scaling variable ^, c.f., Eq. (jSj). Adopting the notation of 
Ref. [32^, the CFFs read for the vector 

Vi^ = -«2-Eif^Wi+(4--PEM + «r-P«y W2-f2-^fi-^5§2-W3 (33) 



and axial-vector 



A. = ^e,...,p7^,- — - (34) 



'A 



4, 



where 1/Q-power suppressed effects have been neglected. 

The CFFs in Eqs. fj33|) and ()34j) are given by a convolution of perturbatively calculable 
coefficient functions C*-^-* and twist-two GPDs via 

nk= ['dxCt\^,x,Q^fi^)Hl{x,7^,A\fi\=.^, for fc = {l,...,5}, (35) 

H/c= J] ^^^xCi+)(e,x,Q^/.2)#^(x,r/,A^/x2)l,=_5, for A; = {1, . . . , 4}, (36) 

i=u,... "^^1 

where n denotes the factorization scale. For each quark species i we have nine GPDs. The two sets 
{HI, . . . , H^} and {HI, ... , HI} are defined by off- forward matrix elements of vector and axial- 
vector light-ray operators. All GPDs of a given set satisfy the same evolution equations, which 
govern the logarithmical dependence on the factorization scale /i. The coefficient functions C^^^ 
have been perturbatively expanded. In LO they read for the even (— ) and odd (-I-) parity sectors 



where Qi is the fractional quark charge. 

The Fourier coefficients c/sj can be calculated from Eqs. ^lU^ . (j7T|) . and analogous ones for the 
squared BH amplitude by summing over the polarization A', where we can employ the common 
projector technique. For a massive spin-1 particle we have (see for instance |15| ) 

et^(A = 0)ei,(A = 0) = S^S„ (38) 



et;.(A = ±l)ei^(A = ±1) = 2 (^~^^^ + Hi^f S^S^ + —e^^p,s 

9 



where A = {+1,0, —1} denotes the magnetic quantum number with respect to the quantization 
direction given by the spin vector 5*^, defined in Eq. (jH)). Obviously, the spin sum of the recoiled 
nucleus is 

J2 ^;,i^')e2u{A') = -9,u + ^^. (39) 

A'=-l A 

As we see, the Fourier coefficients for a spin-1 target quadratically depend on the spin vector 
and, thus, we introduce the following decomposition 

Cn = l^^cl^np + A cos(6') + <tp(<^) MO)} + (^1 " ^A^^ 

X KltpIv') sin(26») + cl^^p cos^O) + <ttp(</') sin^O)} (40) 

for T = {BH,X, DVCS}. An analogous decomposition holds true for the odd harmonics s^. The 
unpolarized and the longitudinally polarized coefficients c/s^^^p, c/s^Lp, and c/s^llp^ respec- 
tively, are independent of ip. The transverse coefficients c/s^ rpp and the transverse-longitudinal 
interference terms may be decomposed with respect to the first harmonics in the azimuthal angle 



c: 



n,Tp(<^) = c^,TP+cos((^) + s^TP-Sin(<^), (41) 
Cn,LTp(¥') = c^,LTP+ COS (y?) +s^^LTP- sin (</?), (42) 

while c/s^rprpp may be written as 

Cn,TTp(<^) = C^,TTPS + cr,TTPACOs(2(^) + S^TTP±sin(2(^). (43) 

Analogous definitions are used for the odd harmonics (just replace c ^ s). Obviously, c^ttps 
does not belong to an independent frequency, rather it can be included in the constant and cos^(6') 
terms of Eq. Indeed, the calculation of the Fourier coefficients establishes the equality: 

cJllp = Kunr> - 2cJttpe, T = {BH, J, DVCS}, (44) 

and an analogous one for rprppj. coefficients. 

We will now present the results for the dominant harmonics c/sf of the interference term. We 
write them as a product of leptonic prefactors C, defined in Tab. ^ with the real/imaginary part 
of 'universal' functions Cj: 

( r^c "1 r ^ 

Cf, for i = {unp,---,TTP-}. (45) 

Since they linearly depend on the nine CFFs and three electromagnetic form factors, we write 
these functions in terms of real valued 9x3 matrices Ma'- 



Cf = (Hf-n^nf-n^^MU G2 \ , for ^ = {unp,---,TTP-}. (46) 











1 4,. J 


H 


Km J 


[ {53m} 
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ric 


rXs 


/-DVCS 


unp 


-SK{2-2y + y^) 


8XKy{2 - y) 


2-2y + y^ 


LP 


-8XKy{2 - y) 


8K{2-2y + y^) 


Al/(2-i/) 


TP+ 




8^^^{2-2y + y') 




TP- 








LLP 


-8K{2 -2y + y') 


8XKy{2 - y) 


2-2y + y^ 


LTP+ 




8X^^^^y{2 - y) 




LTP- 


sVlElMA^2-2y + y^) 


^VlEMAXy{2 - y) 




TTPS 


-8K{2 -2y + y') 


8XKy{2 - y) 


2-2y + y^ 


TTPA 


-8K{2-2y + y') 


8XKy{2 - y) 


2-2y + y^ 


TTP± 


8K{2-2y + y') 


8XKy{2 - y) 


-Az/(2 - y) 



Table 1: Kinematical prefactors that appear in the Fourier coefficients of the interference term 
(gSl) and squared DVCS amphtude (jUj). 

In the case of a longitudinally polarized target the three matrices A^Jnp' -^lp^ ^'^d A^llp 
presented in Appendix IB. II for small values of the momentum transfer^, i.e., — -C M\. 
The squared DVCS amplitude is written in an analogous manner, namely. 



^7"^ = for. = {unp,...,TTP-}, (47) 



where the nine functions 



(4J 



are expressed by 9 x 9 hermitian matrices A^f^*"^. The leptonic factors C^^'^^ are given in Tab.C] 
and the approximation of the matrices for longitudinally polarized target can be found in Appendix 
IB. 21 The analogous BH contributions are listed in Appendix IB. 31 

Let us summarize. For a given harmonic in we have altogether nine possible observables. In 
principle, they can be measured by an appropriate adjustment of the spin vector and Fourier 
analysis with respect to the azimuthal angle ip. The interference term linearly depends on the 
CFFs and is, thus, of special interest. In facilities that have both kinds of leptons it can be 
separated by means of the charge asymmetry. Combined with single spin-flip and unpolarized 

^ The complete expressions for all matrices of the interference and squared DVCS terms as well as 
the complete squared BH term are available as TeX or MATHEMATICA file from the authors via e-mail: 
axel.kirchner@physik.uni-regensburg.de, dieter. mueller@. physik.uni-regensburg.de (dmueller@theorie.physik.uni- 
wuppertal.de). 
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as well as double spin-flip measurements, the sin / cos harmonics can be separated. This would 
provide access to the imaginary and real part of the nine linear combinations Cf and so to the 
nine GPDs. To extract information from the nine GPDs one needs a realistic model for them, 
where the model parameters have to be adjusted by fitting the measured data. Certainly, this 
goal requires a deep understanding of non-perturbative physics and a dedicated facility. 



4 Models for nucleus GPDs 

The main uncertainty in the estimate of observables is the lack of knowledge on GPDs. Since GPDs 
are hybrids of exclusive and inclusive quantities, we may expect that the constraints on models 
arising from the reduction to parton densities in the forward limit or to elastic electromagnetic 
form factors, appearing as the lowest moment of GPDs, ensure the right order of magnitude for 
cross sections and asymmetries. However, one should keep in mind that the predictions for DVCS 
in terms of GPDs are rather complex and that only few quantities, like the sign and the size of 
the single beam spin asymmetry or the size of the unpolarized cross section, are "predictable" in 
a rather model independent manner. We can also rephrase this statement in a positive sentence, 
namely, most of the observables are sensitive to details of the GPD models. One example is for 
instance the charge asymmetry for DVCS on an unpolarized proton target. However, we should 
emphasize again that the interpretation of experimental data requires a careful consideration of 
possible kinematical effects as it has been stressed in j7]. 



4.1 Scalar target 

To start with the simplest case, we want to consider a scalar nucleus target A with charge Ze. In 
this case only one GPD arises at twist-two level and two of them at twist-three level jSHl El El IS] • 
The former one is defined as expectation value of twist-two light-ray operators 

if^/^(x,r/. A, = / ^e-^+-(A(P2)|V^.(-«:n)7+^.(«:n)|A(Pi)) , (49) 

where the gauge link has been omitted. Here the H — component of a four vector is given by the 
contraction with the light-like vector n, e.g., P+ = n-P. Thus, — 1 < x < 1 can be interpreted as a 
momentum fraction with respect to the light-cone component P+, while the skewedness parameter 
7] appears as a scaling variable. Lorentz invariance requires that the nth moment with respect to 
X of this functions is a polynomial in rj of order n + 1. This is ensured by the GPD support, which 
can be implemented in a simple manner by the so-called double distribution (DD) representation 

Hill 

H'/^{x,r],A^,Q'^)= [ dy [ dz x 6{x - y - r]z)h'^^{y, z, , Q^) . (50) 

J-l J-l+\y\ 

Note that, in contrast to the common definition, we included here an extra factor x, which guar- 
antees that the term of order r^""*"^ is not absent in the nth x-moment. Moreover, time reversal 
invariance combined with hermiticity requires that they are even in rj and, thus, the double dis- 
tribution h is even in z. The polynomiality condition is separately ensured for the two regions 
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y > and y < 0. The so-called positivity constraints, mentioned before, are not implemented in 
this representation. 

For A ^ the GPD reduces to the parton and anti-parton distributions 

l^^^^lfx o^i-l H^^Hx,v = o,^' = o,Q') \ .... 

I^/aJI^'^ J 1-/7VA(_^,^ = 0,A2 = 0,Q2)/ ^^^^ 

in a nucleus with the familiar interpretation. Both of them have the support < x < 1 and we 
applied here the standard sign conventions for anti-quarks, so that g*''^ is positive. They satisfy 
the charge sum rule 

= Z, (52) 

and the Baryon number conservation sum rule 

= A. (53) 
In both sum rules the sea quarks drop out and we can read off the number of valence quarks 

/ rfxg"™'/^(x, Q2) = iV„,^, = 2Z + iV , 
Jo 

! dxq'^^^'/^{x, Q2) = N^^^^ = Z + 2N, N = A-Z. (54) 
Jo 



Furthermore, we have the momentum sum rule: 



[\xx W^''ix,Q^)+t'^^ix,Q')]+ f dxxg\x,Q^) = l. 

-^0 i=u4,s 



(55) 



where g^{x, Q^) is the gluon density. 

Now we like to connect the nucleus parton distributions to the proton ones. Isospin symmetry 
breaking effects are rather small and, thus, we employ this symmetry to connect their parton 
content: 

qU/p _ qd/n ^ qU -u/p _ -d/n ^ -u ^d/p _ ^u/n ^ ^d -d/p _ -^/n ^ -d 

qs/p ^ -s/p ^ qs/n ^ -s/n ^ ^ f = g'' ^ Q . 

Considering the nucleus as a system of almost free nucleons, we can immediately express its parton 
distributions by those of the proton 

g«/A _ ^ j^qd ^ qd/A ^ ^^d ^ ^ 

q^'^={Z + N)q\ g^ = {Z + N)g. (56) 
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Additionally, we have to relate the momentum fraction variable x to that of the struck nucleon. 
Since the momentum of each nucleon is the l/Ath part of the nucleus one, the appropriate scaling 
of the momentum fraction gives for instance 

g"/^(a;) = Ae{l - \xA\) [Zq^'^xA) + Nq-^^xA)] (57) 

and analogous for the other parton distribution. Here the loss in phase-space due to the restriction 
X < l/y4 is compensated by the prefactor A. Note that does in reality not vanish at 

X = 1/A, however, this simple prescription is for our purpose justified as long as xb = Ax a is 
in the valence quark region or below. It is easy to convince ourself that all sum rules 



will be obeyed, if they are satisfied for the nucleon. For the ratio of the structure functions F2, 
normalized per nucleon, in deep inelastic scattering one finds for an isoscalar target (A^ = Z) 

— 7 — -7 9 ■i.r, — r = 1 , with xa = Xb A. (58) 

F^{xb) AxBE^Qh'^''i^B) 

Certainly, experimental measurements show that this is only true for xb ~ 0.25 and in other 
regions one finds characteristic deviations from this ratio, which reflect different aspects of the 
binding forces between the nucleons. For the kinematics we are interested in, i.e., xb ~ 0.1, it is 
a small difference of the order of a few percent, which can safely be neglected for our purpose. 

Inspired by the parton picture, we introduce a terminology for the GPDs. However, we remind 
that we deal now with a generalization of distribution amplitudes and the probabilistic interpre- 
tation arises only in the forward limit, because of the optical theorem. We uniquely separate a 
GPD in valence- and sea-like ones 

H'/^{x, T]) = e{x > -|r/|)i^^™'/^(x, 7]) + W^^^/^{x, T]) . (59) 

The sea-like GPD is antisymmetric in x and can be separated in quark and anti-quark ones. 
The latter are obtained from Eq. ()50|) due to the restriction y < and the whole sea might be 
expressed by the anti-quark contribution 

H''^^/^{X, 7]) = ir^^{x, T]) - lt^^{-X, 7]). (60) 

The valence-like part follows by subtracting the sea quark GPD from the whole one: 

i7^™'/^(x, 7]) = e{x >-\7]\) [H''^{X, 7]) - i7'^°^/^(x, 7])] . (61) 

Note that the valence-like GPDs have no symmetry property and that we included the region 
^ 3; < \7]\ in our definitions, although their partonic interpretation belongs more to the 
excitation of a mesonic like state. However, Lorentz invariance requires that the 'inclusive' re- 
gion, which is closer to a probabilistically partonic interpretation, is entirely determined by the 
'exclusive' one - the reverse is generally not true. 

Let us now discuss the form factor aspects of GPDs. We define for any parton species its own 
partonic form factor, which is given by the lowest moment: 

dx H'/^{x, f], A\ Q2) = NiF\A^) . (62) 
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The above given definition induces the equahty of sea- and anti-quark form factors. Moreover, 
Lorentz invariance and current conservation ensure that the r.h.s. is independent on r] and Q^, 
respectively. We normahze the partonic form factors to one 

FXA2 = 0) = 1. (63) 

Thus, Ni might be interpreted as the number of a given quark species i inside the nucleus. 

The above given sum rules for parton densities can now be generalized to the non-forward 
kinematics. The charge sum rule ()52|) turns into the relation between partonic and electromagnetic 
form factors: 

F(A2) = g„iV„^^,F^™i(A2) + g,iV,^^^F'^v.i(A2) , ^(A^ = 0) = Z, (64) 

while the Baryon conserving ones (|53p reads 

^Bar(^2) ^ ^N^^^^p^^^ ^A^) + liV,^^,F'^™' ( A^) , ^^-(A^ = 0) = A , (65) 

The sea-like GPDs are antisymmetric in x by definition and so they can not contribute to these 
sum rules. The first moment of the singlet combination, including the gluonic component. 



1 

1 



i^(A(P2)|T++|A(Pi)) (66) 



n ■ yx, II, ^ , ^ ; -I- -n -■ {x, //, la , '^^^ 

J,=u,d,s 

is given by the expectation value of the H — components of the energy- momentum tensor T^i^. Since 
of current conservation, it is also independent on the scale Q. The r.h.s. is then a polynomial in 
1] of the second order whose coefficients are given by the two gravitational form factors, which 
appear in the matrix element (A(P2)|^^ti/|A(Pi)): 

1 

P 



,2(A(P2)|T++|A(Pi))=Ti(A2)+r/2T^(A2), with T\A^ = 0) = 1 . (67) 



The normalization of T^(A^ = 0) is consistent with the momentum sum rule (j55j) . It has been 
argued in Ref. [3Bj that the second form factor is related to the spatial distribution of strong 
forces inside a nucleus and that its partonic content can be accessed via the GPDs. 

To relate the nucleus GPDs to those of the proton, we rescale, as in the forward case, the 
longitudinal momentum fractions, i.e., their x-dependence and r^-dependence. This is again 
motivated by considering the nucleus as a system of almost free nucleons, which have approx- 
imately the same momenta Pi/A = (1 — ri)P+/2A. At leading order in the interaction 
with the two photons take then place on one of them, which has the outgoing momentum 
Pi/A + A = (l-r/ + 2At])P+/2A + A^. Thus, instead of the GPD definition ^ we have 
to consider the matrix element, see Fig. 

e'^P.^N{P,/A + A)m-Knh^ij,{Kn)\N{P,/A)) . (68) 

To make contact with the nucleon GPD, we refer in this matrix element to the light-cone com- 
ponent of the struck nucleon P^ = [1 — 77 + ArjjP^/A. This implies the rescaling prescription for 
the light-cone momentum fraction and skewedness parameter 

xn = \ — —Ax and t]n = ^ = \ — —Ar] . (69) 

1 — ?7 P" 1 — 7] 
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Ai dhk 

(a) (b) (c) (d) 



Figure 3: Different contributions to a nucleus GPD in which the hght-ray operator couples to a 
nucleon (a), meson exchange in the t-channel (b), mesonic contact term (c), and virtual nucleon 
pair (d). 



From the second equality one easily recovers the scaling law for the kinematical variable, namely, 
xb = Ax A- However, this picture can not be quite realistic, since transferring the whole transversal 
momentum A"*- to the struck nucleon, should, compared to a binding energy of about 10 MeV 
per nucleon, induce target dissociation. On the other hand, nuclei form factors are non-zero even 
for rather large momentum transfer, which shows that there is a certain probability that the 
transversal momentum transfer is almost equally distributed among the constituents^. Assuming 
that the A^-dependences factorize, we can take care of this bound state effect by replacing the 
partonic proton form factors with the nucleus ones: 



i7'^-'/^(x,r/, A^) = F'^vai/A(^2 
A^) = F^<^=^/A(a2) 
i/^/^(x,r/. A') = Ff/^(A2 



dXN 






dx 


0{\xn\ 


dXN 




dx 


0{\xn\ 


dXjy 




dx 


9{\xn\ 



< 1) [Zff'^™' + OT"-'] {xn,Vn) 



(70) 



dx 



N 



dx 



9{\xr,\ < l)AH3{xM,m) 



Here we introduced the reduced proton GPDs by W{x, rf) = H\x, r], A^ = 0). The support of the 
proton GPDs induces the restriction A\x\ < {1 — r])/{l — i^n), which of course is only an artifact 
of our approximation. It is easy to check that by construction all sum rules remain valid. For an 
isoscalar target isospin symmetry requires 



ri, A^) = H'^/^{x, ri, A^ 



(71) 



and, thus, also the partonic u and d form factors are equal. 

Obviously, what we presented here is the simplest version of a convolution model for nuclei 
GPDs, which certainly can be improved by incorporation of binding effects. Hopefully, they 
might provide only small corrections to the scaling law (f7n|) and should explain the distribution 
of transversal momentum transfer. Nevertheless, in the case that such effects are responsible for 

^Certainly, it remains an experimental problem to ensure that the measured Compton scattering process is 
coherent and a theoretical task to derive the predictions for target dissociation. 
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symmetry breaking, e.g., the breaking of spherical symmetry imphes a non- vanishing quadrupole 
moment for deuteron, they might become a leading contribution to GPDs that decouples in the 
forward case. Thus, it is worth to have a closer look on them. 

First we consider the possibility that the light-ray operator insertion couples to meson exchange 
currents in the t-channel, as it is shown in Fig. [Hb- The corresponding matrix element reads 



fj 

/TT 



-Kn)^+iJi{Kn)\p{PjA)n{P^/A)) 



t— channel 



oc / du(f){u) 



-1 



dx 



dn 
— ( 
27r 



{n{[u - 1]A/2)|0(k, -i^M[u + l]A/2)) 



.^4 



Here the probability amplitude (f){u) with |m| < 1 is a function of the light-cone momentum fraction 
u, defined by P^^ = nA_|_ = urjP^. The momentum fraction x"^^ and skewedness parameter 77®^, 
entering in the GPD, are 



x^ 



X 



urj 



and rf 



PI 



1 

u 



The support property of the resulting GPD 



i/*-~'(x,r/) oc f- / dusign{u7])(t){u) 

\m J-i 



X 



1 n-\y\ 

dy dz 5{x — rj[uy + z])F^ 
1 J~l+\y\ 



t— cannel 



(72) 



(y, ^) (73) 



can be deduced from the DD representation. It only contributes in the 'exclusive' region, since 
for the 'inclusive' region we find \x/ri\ > 1 > \uy + z\. Time- reversal invariance combined with 
hermiticity requires that (p and ^^-c^nnoi^^^ symmetric in u and z, respectively. Consequently, 

the exchange GPD is also symmetric in x and the polynomiality condition is satisfied in general. 
For the lowest moment we have the normalization: 



dxH 



t— cannel / 



{x,ri) (X / du\u\(t){u) / dyy 



^-\y\ 



dzF^ 



t— canncl 



(74) 



Remembering that the integral over z gives the parton densities of the pion, i.e., qn{y)/y, we 
obviously, establish a link between 'exclusive' t-channel contributions and the parton picture. 
Again the total sea quark contribution q^^{y) is antisymmetric in y and, thus, it drops out in 
the sum rules fl64l65p and also for higher moments (which are even). Such an J = 1 isospin 
exchange contribution is 'filtered' out for an isoscalar target, however, it might be important for 
other targets or dissociation processes. Note that for an isoscalar target, resonance exchange 
contributions appear, like the famous pirj exchange current. We should add that in such a case 
also contact terms, depicted in Fig. IHt, have to be included, which for instance are given by 



dn 
2^' 



(iV(PiM + A/2)7r(A/2)|V^,(-Kn)7+7/>,(/s:n)|iV(Pi/A)) 



(75) 



Moreover, Lorentz invariance requires that virtual nucleus anti-nucleus states contribute to the 
nucleon GPD. An example is depicted in Fig. Ell and its matrix element reads 



fj 14^ 

ZTT 



^4 



(76) 
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It is far beyond the scope of this paper to estabhsh a closer hnk between the effective forces of 
nuclei and the fundamental degree of freedom in QCD, however, our discussion shows that GPDs 
are an appropriate tool for this task. 

Based on a 'popular' model for the reduced proton GPDs, we present now our model for a 
scalar nucleus. To implement the support properties we deal with the DD representation ^1 12] at 
a given input scale Q = Qo- 

/I r^-\y\ 
dy I dz 6{x — y — riz)q\y)iT{\y\, z;b^) + sign{ri)D\x/ri) (77) 
■ 1 J-i+\y\ 

and adopt Radyushkin's proposal for the factorization of the DD into a forward parton density 
q^{y) given at the input scale Qq and a profile function 



which is normalized to 1: 



dzTT{\y\,z;b) = 1. (79) 



-n 



Here the free parameter b models the skewedness effect that arises from the x-shape of GPDs for a 
given value of t]. In the limit b ^ oo the GPDs are independent of r] and are simply given by q'^{x) 
itself. For the unpolarized parton densities of the proton we will take the MRS A' parameterization 
|46j . given at the input scale Qq = 4: GeV^, with u = d = s/2. For the electromagnetic form factors 
we adopt here for simplicity the following parameterization 

z 

where the slope 7^2/3 j^sea ^^^e sea quarks is defined in such a way that B^'^^ only weakly depends 
on A. This is motivated by the fact that the nucleus density is rather independent on A and, 
thus, we expect that the radius of the sea quark distribution should scale with A^^^. For the 
electromagnetic form factor we adopt the parameterization of Ref. |47j: 

1 



F^(A2) = Z (^1 - ^c^A^^ exp I^A'^'I 



with b = 2.04 X 5.07 GeV"^ and c = 1.07^^/3 ^ 5.07 GeV"^ 

There are a few comments in order. First the factorization in A^ for medium values has 
no theoretical or experimental justifications. In comparison with experimental data, which are 
presently measured at a mean value of < — A^ >~ 0.1... 0.3 GeV^, the A^-dependences effectively 
enters only the normalization of the CFFs via the slope parameters. Second remark concerns 
the polynomiality property of the moments. In comparison to Eq. (j50|l we choose for simplicity 
another DD representation in Eq. (|77jl . The prefactor x is now neglected and the highest possible 
term r/"+^ for the nth odd x-moment is restored by the so-called D-term. For detailed discussions 
on this subject see Ref. |48| IE] . This D-term is antisymmetric in x and can not spoil the sum rule 



18 



however, it enter the sum rule ^7\i and is discussed below. In our language we count it as 
a sea-quark contribution. The third remark concerns the factorization of the DDs ansatz, that 
is inspired by an intuitive picture and provide in general an enhancement at the point x = ±1]. 
Finally, we should comment on the meson exchange contributions. For a scalar target we will 
not explicitly take them into account. They effectively enter here the parameterization of the 
form factors. Moreover, they have do vanish in the forward limit and, therefore, they should be 
suppressed for valence like contributions at small A^, while for sea-like contributions they can be 
included in the antisymmetric D-teim. 

Let us finish this section by expressing the nucleus CFF Ti.^ by the proton ones in the twist-two 
and twist-three sector. From the ansatz (ffn|) and the convolution^ ()H5|I . one easily establish the 
following scaling law for the twist-two CFF of a isoscalar nucleus 



2^ ' '\Ql Q 



(82) 



Here ^{^n) denotes the contribution that appears in the decomposition of the proton CFF 
Ti = 7i"™' +7i'^™i +'H^^^ at = 0. Note that, as long we rely on the ansatz (fTOI) . this scaling law 
is valid beyond the LO approximation, since the hard-scattering coefficients have the following 
functional dependence l/^C{x/^). It is interesting to note that in Ref. it was argued that the 
D-term contribution to the Compton form factor, which is independent of xb, scales with A^^^. 
Its contribution to the second gravitational form factor is 

T^(A^) = i|^ E l\xxW/\x,v,A') = U'd\A'). (83) 

i=u,d,s,g ^ 

To study the D-term with respect to A and in relation to the sea quark contribution per nucleon, 
we have rescaled by A"^. Based on the estimate in the same reference, one finds that ^^(A^ = 0) 
slightly increases with A > 8: 

d^{A^ = 0) = rf«/^(A2 = 0) + d^/^{A^ = 0) ^ -0.2^1/3 ^1 ^ 3.8/^2/3) , (84) 
rfQ/^(A2 = 0) ^ -4.0 at /i ^ 0.6 GeV, (85) 

where the second result for the nucleon D-term in the quark sector has been predicted by a model 
calculation within the chiral soliton model ED]- As it will be discussed in Section 15.11 a 
present measurement of the charge asymmetry on the proton target does not allow for a definite 
conclusion of the D-term contribution. Unfortunately, we might conclude that for nucleus target, 
e.g., ~ —1.3 for A = 200, this term will not induce a significant contribution, too. Certainly, 
it is an important task to confront this expectation with experimental measurements. 

Employing the equation of motion, the twist-tree GPDs are decomposed in the so called 
Wandzura-Wilczek (WW) term, entirely expressed by the twist-two GPD H, and the quark- 
gluon-quark GPD H^'^^, carrying new dynamical information. The two twist-three GPDs enter 
in the DVCS amplitude always in the same linear combination [H]: 

^cff _ ^cfT-WW ^ ^qGq^ (gg) 



^Remove the index k for a scalar target 
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where 



Using the ansatz ()7()|1 . one easily derives the WW piece in terms of the proton GPDs H^. For 
instance, the total sea quark contribution reads in terms of the SU(2) symmetric antiquark GPDs 



2e 



1 + C 1 + ^ ^N + x ^N-x-tO 

t=u,d,s 

x2^^'^{x,iN) -li'""{-x,iN)] , (88) 

where = ^^/(l + ^ — ^0 is to be considered as a function of C,. The analogous formula in 
the valence quark sector follows by an appropriate replacement of form factor and GPDs together 
with the corresponding charge factors. 

4.2 GPDs of spin-1/2 nucleus target 

Analogous models for proton GPDs have been used in Ref. [7], for a detailed discussion of GPDs 
see Ref. j^, and they are consistent with all experimental data available at present in the LO 
analysis. It has been stressed in Ref. [^2] that this not necessarily true in NLO. The problem 
arises from the employed model of the flavour singlet GPDs, especially from the gluonic ones. We 
will skip this issue here and refer to the discussion in Ref. [7]. Due to the spin structure of the 
nucleon, one has to model four different sets of GPDs: 

F=[h,E,H,e], (89) 

corresponding to helicity conserved and non-conserved form factors with even and odd parity. 
Note that for an unpolarized target the dominant contribution arises from the helicity conserved 
and parity even GPDs H^. They satisfy the constraints 

\imW{x,r],A'^)=q\x) and / dx H\x , r] , A^) = Fi{A^) (90) 

For longitudinally polarized target both H and H are important. The latter ones are related to 
the polarized sea quark content and axial-vector form factors: 

\imW{x,r],A^) =Aq'{x) and [ dx H\x , r] , A^) = G\{A^) , (91) 

A^O J_i 

respectively. The helicity non-conserved GPDs E and E, which decouple in the forward limit, 
obey the sum rules 



j dxE'{x,r],A^) = Fl{A^) and j dx E'{x,r], A^) = Gi{A^) . (92) 
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Perhaps they can be accessed in transversally polarized target experiments. An essential obser- 
vation was that the ansatz (f77|) for the sea-quarks, measured in the small xb region, has to be 
suppressed by choosing a b^"^^ oo and a large slope B^'^^ ~ 9 GeV^. This suppression then 
ensures also the correct size of the single beam spin asymmetry measured in fixed target exper- 
iments. For instance, the mean value at HERMES is xb ~ 0.1 and, thus, the value of ^ ~ 0.05 
appearing in the argument of GPDs is rather small. As we said, the contribution of sea quarks 
remains relatively small with respect to the valence like ones. If this observation is implemented 
in the correct way within our GPD model remains an open issue at present. 

The helicity conserved GPDs H^^^ and H^^^ of a nucleus with spin-1/2 can be obtained from 
those of a proton by means of Eq. (ffn|) . Unfortunately, not too much is known about the structure 
of E'/^ and E'^^, beside the pion-pole contribution of For the parametrization of the twist- 

three GPDs see Refs. 

4.3 GPDs of spin-one target 

In this section we list the properties that are known about the spin-one GPDs and have been 
already given in Ref . • For the matrix element of the twist-two operators we use the analogous 
kinematical decomposition as for the CFFs in Eqs. (|33|) and (|34p . Combining hermiticity and time 
reversal invariance tells us that the GPDs are real valued functions, which respect the symmetry 
properties: 

Hkix,!]) = Hk{x,-r]) for = {1,2,3,5}, Hi{x,r]) = -Hi{x, -r]), 

Hk{x, f]) = Hk{x, -f]) for k = {1, 2, 4}, Hs{x, r]) = -Hs{x, -t]). (93) 

Note that Hi{x,rj) and H-^{x,ri) are antisymmetric with respect to 1], and, consequently, their 
moments are odd polynomials in 7]. 

From the definition of GPDs as Fourier transforms of light-ray operators it follows that their 
lowest moment is given by the form factors appearing in the vector or axial- vector current, re- 
spectively. In the former case they are related to the electromagnetic form factors 

Y,Q^ j'^dxH'J\x,v,A\Q') = Gk{A') for A; = {1, 2, 3}, (94) 
j dx H'J^{x, T], A\ =0 for A; = {4, 5}, (95) 

or they have to vanish. For the deuteron the electromagnetic form factors are known from exper- 
imental measurements and their parameterization is given in Appendix |X] [SH EHl- For a target 
with charge Ze we choose the normalization 

Gi(A2 = 0) = Z, G2(A2 = 0) = ZfiA , GUA^ = 0) = ZifiA + Qa - 1), (96) 

where Zha is the magnetic moment and ZQa is the electrical quadrupole moment. The latter 
is induced by non-central nuclear forces. Its non-vanishing value, measured for the deuteron, has 
been considered as an evidence for the role of pions in nuclear physics. More precisely, the one 
pion exchange provides at larger distances the dominant contribution to the potential, which is of 
the Yukawa type, and induces a D-wave admixture to the S-wave function. 
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The sum rules for parity odd GPDs read 



l\xHl{x,v,A',Q') = Gl{A') forA; = {l,2}, (97) 
j dx Hi{x,r], A'^, Q'^) = for k= {3,4}. (98) 

The two axial form factors G^(A^) with k = {1,2} are defined by the matrix element of the 
current j^'* = ipi{0)'y^'y^ipi{^) with flavour i: 

jr = fe,.5..pG, ( A=) - '^>.^P..^2-^^y>.AP.5^.-^ g^(A.) , ,99) 

In principle they, i.e., certain linear combinations with respect to the flavour number, can be 
measured due to weak interaction, however, to our best knowledge this has not be done yet. 

First we consider the forward limit in which H2,H^, H4 as well as H2,H^, H4 decouple from 
the Compton amplitude and, thus, are not measurable in deep-inelastic scattering. This does not 
mean that the functions in question vanish by themselves. The remaining three functions are 
expressed in terms of parton densities 

q'/^{x) = Hl^^ix, r] = 0,A^ = 0) = l- {q+\x) + q-\x) + q%x)} , 

o 

Sq'/^ix) = H5{x,ri = 0, = 0) = g°(x) - ^ {g+'(x) + q-\x)} , (100) 
Ag'/^(x) = Hi^^{x,r] = 0, A^ = 0) = g+^(x) - gf^(a;), 

where q^ = q^ + q^. Here qf{x) is the probability to find a (anti) quark with momentum fraction 
X > {x < 0) and positive helicity in the target of helicity A. Note that these definitions contain 
both quark {x > 0) and antiquark [x < 0) contributions with the following sign conventions: 

q{x) = —q{—x) , 5q{x) = —6q{—x) , Ag(x) = Ag(— x), for x > 0. (101) 

The combination of quark distributions in H5 enter the structure function bi, measurable in deeply 
inelastic scattering on a polarized spin-one target. The sum rule then induces 

1 .1 

dx6q{x) = =^ / dx6q'"'\x) = 0. (102) 
-1 Jo 

To obtain the second sum rule for valence quarks we employed the antisymmetry property of the 
sea quarks. This relation then converts into a sum rule for 61 



[ dxbiix) = 2 / 



{x)dx, (103) 



which vanishes for an unpolarized quark sea . Considering the deuteron as composed of almost 
free nucleons induces a vanishing ratio [28j. The first preliminary measurement from the 

HERMES collaboration j^Tj indicates a tensor asymmetry, A^z = — 26i/3Fi, that strongly depends 
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on xb- It is compatible with zero in the valence quark region, but significantly positive and negative 
for large and small xb, respectively. Since we are in the following interested in the valence quark 
region, we set Sq^{x) = 0. Consequently, this simplifies the spin content of the unpolarized quark 
distributions 

gO(a;) = l{g+i + g-i}(x) ^ qix) = ^ {q^' + q~'} (x) . (104) 

Let us now discuss the modelling of GPDs. In this paper we only consider numerical estimates 
for unpolarized or longitudinally polarized target, in which H2, -^4, H2, H^, and H4 are relatively 
suppressed by ^ ~ (xb/2A) or r = A^/4y4^M^ in the DVCS amplitude and will be neglected. To 
satisfy the sum rule ()95|) for H^{x, 77) we set it to zero. This choice is motivated by the knowledge of 
parton densities. Alternatively, we can also take an antisymmetric function in x, e.g., proportional 
to Hi{x,ri) — Hii^—x,!]). It is worth noting that both GPDs can be related to the proton GPDs 
E and E, respectively. They are accessible in experiments with transversally polarized targets, 
which, however, is beyond the scope of this paper. 

For the GPDs H^^ and Hl^^ we employ the scaling relation ()70|) to connect them to the proton 
GPDs W and H^, cf. Eqs. fj56|) and ()57|) . For the latter we will neglect the D-wave admixture, 
which provides in the forward limit for the deuteron a few percent effect. To implement the 
support properties of GPDs as well as the reduction to the parton densities in a simple manner 
we again employ the DD representation 

Hi} (G\] , f^~\y\ , ^ f q\y)n{\y\,z;b\) 



, . , (A') / dy dz 6(x - y - r]z) < -^-^^'-"'i k (105) 

HI} [G\}^ 'J-iJ-iMy\ ' '\Aq^{y)7r{\y\,z;b\)S ' ^ ' 

where a possible D-term is neglected. The partonic form factors are normalized to one and for an 
isoscalar target we choose 

^n.al(^2^ = ^^^'(A^) = -^(A^) . (106) 

Analogous to the scalar case, cf. Eq. this choice ensures that the sum rules are satisfied, 
where again the sea quarks do not contribute. Inspired by the counting rules, we take for the 
deuteron a partonic sea quark form factor with a rather simple parameterization 

Gr(A2) =(1 + '^Q'^^] (1 - , (107) 



3 /V rn: 



sea 



where the free parameters 7«sea = Q/B^^^, expressed by the slope B^'^^ of the 'charge' form factor, 
and the 'quardropol' moment for sea quarks are specified below. Although theoretical consider- 
ations predict also a rather complex shape for the axial form factor, see Ref. [SH] and references 
therein, we will use the simple parameterization 

G\{A^) = (^-^^^ for z = K^\d™'}, (108) 
Gi(A2) = ( 1 - ^A' I for i = rf^"", 



23 



for the partonic form factors in the small — region. The normalization G\{A'^ = 0) = 1 ensures 
the correct reduction to the forward limit. 

As mentioned above, the quadrupole moment arises from a pure bound state effect, successfully 
explained by the one pion exchange, and essentially given by the matrix element sandwiched 
between S- and D-wave state of the deuteron. As stated before, exchange and contact 1 = 1 
contributions drop out and so the dominant contribution should arise from the overlap 

Hf- [ ^e^^''^^diP2Mm-'^n)^+U'^n)\d{P,,%)) (109) 

Certainly, we have to include the interaction between both nucleons. For the integrated GPD, 
i.e., the form factor G3, one can write the form factor as a product of the isoscalar nucleon form 
factor and the so-called body form factor. The latter mainly arises due to the overlap of D and S 
waves and alters the dependence of the isoscalar form factor. An extension of this approach 
to GPDs has been proposed in j^S], where is essentially given as convolution of the isoscalar 
GPD, modelled as = {H^ + H'^)/2. On the other hand, is probing the binding force of the 
deuteron, and as we know from deep inelastic scattering this effect rather depends on the value 
of x-Q- For GPDs the situation is rather complex and at present we have no calculation or model 
available that allows us to have a deeper insight in this problem. Thus, we consider two extreme 
cases and discuss their consequences in the next section: 

• We assume that arises from the binding forces between 'partons' that carry either rather 
large or small momentum fractions. Thus, one might expect that such effects only slightly 
influence DVCS observables at present fixed target experiments. Therefore, we set for mod- 
erate values of x-q the OFF to zero: 

7^3 = for 0.1<xb<0.3. (110) 

• In the convolution model Hi and are both essentially determined by the isoscalar GPD 
H^^°. So we equate the reduced GPDs Hi and H^: 

H^{x,r,) = Hi{x,'n) . (Ill) 

5 Estimates for observables 

In Section 15.11 and 15.21 we give analytical and numerical estimates for the size of the unpolarized 
cross section, the beam spin asjnnmetry 

^^^^'^^-daT(0) + ciai(<^)' ^^^^^ 
the longitudinally polarized target spin asymmetry 

da^{<t>) - da\4>) ..... 
^"^^'^) = dan0) + rfa^(0)' ^'''^ 
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as well as for the charge asymmetry of the unpolarized cross section 



^C(0) 



da+{(f)) -da~{^) 



(114) 



da+{(f)) + da- {(f)) ' 



Moreover, for a spin-1 target we consider also the beam spin asymmetry 



^L±(0) 



da^^{(P) + da^^{(f)) - da^^{(f)) - da^^{(f)) 
dan{(j)) + rfat^(0) + dal^{(j)) + dai^{(t)) 



(115) 



and the tensor polarization 



da- 



zz 



da^z = da^ + da^ — 2da"^ , da^^^p = - {du^ + da^ + da^^ 



(116) 



3 da. 



unp 



where i[, Jj-, and =^ denote the magnetic quantum numbers A = {+1, —1,0} for a longitudinally 
polarized target. We give estimates for HERMES and (upgraded) JLAB kinematics for a lepton 
beam of E = 27.6 GeV and E = 6 (12) GeV, respectively. 

The single spin and charge asymmetries are dominated by the first harmonics of the interference 
term (f?fj) . which arise at twist-two level, while for the tensor polarization the squared BH term 
dominates in the charge even sector. All observables are contaminated by 0(1/ Q) suppressed 
effects, induced by other harmonics in both the interference and squared DVCS term or from 
pure kinematical effects: higher harmonics of the squared BH term ()25p in the denominator and 
additional dependence due to the BH propagators ViV2{(p)- We found that higher twist-three 
harmonics are small in the WW approximation [7]. For the charge asymmetry also a 'constant' 
contribution arises in the interference term, which is completely determined by the twist-two 
GPDs. It can be relatively large with respect to the dominant cos((/>) term, since in contrast 
to higher harmonics it does not vanish at the kinematical boundaries. Certainly, both the size 
of asymmetries and their power suppressed contamination depend on the ratio of DVCS to BH 
amplitude, which is estimated to be'' 



For HERMES kinematics with jA^] < 0.3 GeV^ xb ~ 0.1, and ^ 2 GeV^ i.e., y ~ 0.5, the 
ratio is T'Dvcs^'^bh ^ 0,^A/Z. For JLAB the value of y is typically larger. This results into 
a stronger suppression of the DVCS amplitude, even for a smaller photon virtuality. We realize 
that for HERMES kinematics the ratio of interference to squared BH term is not necessary small, 
which is obviously consistent with the measurement of sizable beam spin and charge asymmetries. 
However, we have maybe overestimated the ratio ()117p . If the sea quark contribution is not yet 
dominant in this fixed target kinematics, we might have an additional suppression by a factor 
y/xB- Indeed, the measurement of the charge asymmetry for a proton target shows that the real 
part of the interference term is suppressed by a factor of about ten. The unpolarized squared BH 
amplitude, which can be taken into account exactly, is dominated by its 'constant' term, while 
higher harmonics are suppressed by the factor K ~ a/— (1 — y)^"^ / Q^- 

^This is based on the assumption that the CFF scales hke A^/xb as it foUows from Eq. (|70|l together with the 
small Xb behavior of the sea quarks. In the valence quark region we expect even a larger suppression. 




(117) 
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Since we expect, in the kinematics we are interested in, that the unintegrated asymmetries 
()112l 11131 I114|) are dominated by the zero harmonic of the squared BH and first ones of the 
interference term, the dependence of the BH propagators will almost cancel each other. Thus, 
their Fourier series read 

Alu(0) = sin((/.;)4\; + ■ ■ ■ , Aul(</>) = M<f>'^)A'il + ■ ■ ■ , 

AM = aP + cos(0;) Ag^ + ■ ■ ■ ' (118) 

where the Fourier coefficients are proportional to certain hnar combinations of twist-two CFFs, 
contaminated by 1/Q^ power suppressed contributions. For the charge asymmetry Ac{(f)) we 
included also the 'constant' twist-three contribution, discussed above. To coincide with the def- 
initions, used by the HERMES and CLAS collaborations, we wrote the Fourier expansion with 
respect to the azimuthal angle of the real photon = vr — for a frame in which the 2;-axis points 

in the direction of the virtual photon momentum. This convention affects only the sign of ^4^''. 
The asymmetries of (p'^ integrated cross sections, e.g., 

^ _ 2/,^-rf0; sin(0;) {da^-da^) _ 2 ^ d<P'^ sin^(0;)/nP2(0;) ^^g) 
t d<P'^ (da^ + dal) ^ STd^'.y^^'^^W 

is influenced by the angular dependence of the BH propagators and is generally smaller as the 
lowest harmonic of the Fourier expansion ()118|1 . Moreover, the dependence on the kinematical 
variables can be altered. Note also that for xb, A^, or integrated asymmetries, the mean value 
for different parts of the leptoproduction cross section can be quite different. 



5.1 Analytical approximation in the valence quark region 

The lengthy expressions for the Fourier coefficients make it rather hard to interpret future mea- 
surements in terms of GPDs. Thus, it would be instructive to have an approximation at hand 
that serves as a guide for adjusting model parameters by using efficient numerical codes. This 
can be done for the observables and kinematics, we discussed above. For instance, the beam spin 
asymmetry ()112j) can be simplified to 

^Lu(0)~±-4Fsi^(</') with l+f'"^, (120) 
y c^^n^ I - for e+ 

where we neglected possible contamination of the squared DVCS term. Moreover, for — Amin <^ 
— A^ ^ M\ and x-q = Axa ^ 0.3 the Fourier coefficients can be drastically simplified due to 
a rough approximation of kinematical factors, i.e., we expand the Fourier coefficients Cp^^p and 
"^f,unp (sse Appendix ini) to zero order in xa and r. However, we have to pay special attention on 
terms that contain and 7^3, since G3(A = 0) is enhanced by one order of magnitude compared 
to the normalization of the other two form factors. Thus, we take also tG^ and rTi^ into account 
and obtain for the beam spin asymmetry: 



2-2y + y^ 



^Lu(0) ~ ± . J , ^ ; (121) 



2^171:1 + (d - 2TG^){rLi - ^tHs) + irG^n 



"^"^ 2Gi + (G, - 2rG,Y ""^'^^ " 
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To demonstrate that these considerations are more on a rough quantitative level, we mention that 
the kinematical prefactor 1 < (2 — y)/(2 — 2y + y"^) < (1 + v^)/2 should for consistency also be 
neglected, which can induce a reduction of 21%. In the following we take it into account. One 
realizes that the single beam spin asymmetry is proportional to a linear combination of nucleus 
GPDs multiplied with xa = x-q/A. Since for large A the form factors scale with Z and the 
CFFs with A?, this asymmetry does not scale with A, however, besides bound state effects it 
only depends on the ratio A/Z. It is rather surprising that bound state effects enter already the 
interference term for the unpolarized target not only due to a possible modification of the scaling 
relation for Tii, but also directly by the appearance of and 7i^. 

In the case that the beam spin asymmetry is measured on an incomplete unpolarized target, 
where only over the polarization states with A = ibl is summed, the prediction differs from the 
asymmetry ()12ip : 



~ ±— BH _VBH Sin(0) With _ 



and reads in our kinematical approximation 



2-2y + y2 

^ 2Gl(7^l-|7^5) + 2(Gl-2rG3)(^l- 2x7^3-1^5) . , 
2GI + 2(G, - 2.^3)^ ""^'^^ • 



^L±(0) ~ ± — , _ ; — (122) 



The approximation of the unpolarized charge asymmetry ()114|) 



Ac(0)^^^'--+i-^^"^^'^^--- (123) 



y c, 



OjUnp 



is obtained in the analogous manner. The dominant twist-two harmonic is given as the real part 
of the same linear combination of GPDs as in the beam spin asymmetry ()121|) : 



Ac{<P) (124) 

y 

^ 2GiHi + (Gi - 2rG3){ni - 2Tn,) + IrG^H^ ^ 

2GI + {G, - 2rGsr '"^'^^ ^ ' ' ' ' 

Also this asymmetry is nearly independent on A. As noted before, the constant term only depends 
on twist-two GPDs and does not vanish at the kinematical boundaries. In our approximation it is 

,125) 

For a scalar target the exact relation does not depend on GPDs jO], while for a spin-1/2 target 
only a xa suppressed dependence appears jHSlEI- For a spin-1 target the analogous relation is 
unknown at present. 
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The single spin asymmetry ()113|) for longitudinally polarized target 

AuL(0)^±^4^sin(0) + --- with (126) 
y Co,unp I - for 



is mainly governed by Tii. 

A (^^ ^ ^^^^(l-y) 3(Gi - rG3)?^i + ^^G2in^ - rH, - H,) . 

^ul(0) ~ ± S>m 2Gf + (g, - 2rG3)^ ^'''^ 

Since the polarized quark contribution is relatively small compared to the unpolarized ones and the 
latter increase faster with growing 1/xb (in the small xb region), we included also a xa suppressed 
term that is proportional to Tii — tHs — H5. 

For a spin-1 target another quite interesting observable is the tensor polarization asymmetry 
pifij) . Their measurement requires a longitudinally polarized target. For an unpolarized beam 
this asymmetry is dominated by the squared BH contribution, which provides a dominant constant 
term that is proportional to tG^: 



A,, = — - — ^ ± 128 

2G2 + (Gi - 2rG3)' y 



x3Gi 



cos(0) X . . . + ^ sin(0) X 

2 - 2?/ + 



+ ■ 



Here the ellipsises, proportional to the first harmonics, arise from a Taylor expansion with respect 
to — A^/ and stand for the real and imaginary part of a certain linear combination of CFFs: 



The third term comes from the denominator, while the cos(0) term that stems from the BH Fourier 
coefficients c^^^^ and cf do not contribute in our approximation. Consequently, to extract the 
interference term the squared BH contribution has to be subtracted. Alternatively, we can form 
the tensor asymmetry from the charge odd part of the cross section 

, , dal — dal a;A cf.unp ~ '^,llp /,n , /ion\ 

Aczz{(t)) = e+ I e- BH COs(0) + ■ ■ ■ , (130) 

where the ellipsis include a constant term and higher harmonics. Another possibility is to use a 
polarized beam and to form the asymmetry 

A (.^\ - ^^^'^^ = ~ da,,{\ = -1) 

^^^^^'^^ 3rfa,,p(A = +1) + 3rfa,,p(A = -1) ^'^'^ 



Xa -sLnp - ■sf LLP . / ,N , .^1 / + for e" 
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which is dominated by the sin(0) harmonics of the interference term only. These two asymmetries 
allow to explore the real and imaginary part of the following linear combination of GPDs: 



11 -(th (132, 



Aczz}' ' V ^ " i-l/yj icos(0) 

53m \ Gi (tHs + n,) + rGs {H, - 2r7^3 - i^s) 



3fJe J 2G? + (G'i-2rG 



3 



\2 



It is now instructive to compare these results with those for a spin-1/2 or -0 target and to 
explore its dependence on GPDs, which are of course model dependent. The analogous predictions 
simply follow by setting G3, TC3, and to zero and by the replacements: 




for spin- 



1/2 




First we like to remind that the sum rules, e.g., Eq. (jMj) . suggest that the A^— dependence of the 
valence-like GPDs is given by Gi and Fi (or the scalar form factor F), respectively. Furthermore, 
the analyses of the HI, HERMES and CLAS DVCS data for a proton target in terms of the 
oversimplified model, given in Section El at LO indicate that the unpolarized sea quark contribu- 
tion is compared to the forward case additionally suppressed^. We suppose that this qualitative 
property holds also true for the nucleus GPDs. Therefore, we neglect in the valence region the 
sea quark contribution and model the remaining GPDs as a product of form factor and valence 
quark distributions^. In LO of perturbation theory we find for nuclei targets 



expressed in terms of the quark densities of the proton. The analogous formula holds true for the 
spin-1/2 case and also for the Hi contribution of spin-1 targets. 

The beam spin asymmetry for the positron scattering off a proton target is 



(134) 



The ratio of beam spin asymmetries for a spin-0 or -1/2 nucleus with Z ^ to the proton is 



^Lu(0) Qlq^-i^N) + Qy-^ii 



N 



for J ={0,1/2}. (135) 

c 



®We remind that £^ ^ xb / {"2 — xb) and already for ^ ^ 0.05 one finds in deep inelastic scattering an almost equal 
contribution of valence and sea quarks. 

^Here we neglect any skewedness effect. For the DD model, introduced in Section^ skewedness effects would 
provide an enhancement of the CFFs. Strictly spoken, it is not known what caused the observed suppression of sea 
quark contributions and thus the simple ansatz for the DD's together with a rather large slope parameter might 
be questionable. 
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For a spin-1 nucleus we would have the same prediction, if 7i™^ = G?JGihC^^ holds true. However, 
if 7i™^ is negligible in this kinematics, it follows from Eq. ()121|) an additional A^-dependence, which 
is crucial: 



Alu(0) 2 + (1 - 2rG,/G,Y Qlq^^-{^N) + Q^-^i^N) 
At^{<P) 3 - 3rG3/Gi (Ql + Qj) {g"-(eiv) + q'-H^n)} 



for J 



Alu(0) 2 + 2(1 - rG,/G,y QIq^^H^n) + Q^^HO 



(136) 



Taking xb ~ 0.1, we expect from the parameterization of parton densities that (ivai/^vai ~ 0.5. 
For A2 = -0.2 GeV^ ^ 2.5 GeV^ and E = 27.5 GeV we estimate the beam spin asymmetry 
for a proton target and the ratios ()135|) and ()136|) for nuclei [Z ~ N) to be: 

^ -0.29. {f } fo. {f >} . (137, 

Our estimate for the sin(</))-weighted asymmetry (jll9|) is Ai^y ~ —0.23, while the ratios of beam 
spin asymmetries for different targets are unchanged. So we roughly expect an enhancement of 
the beam spin asymmetry for nuclei with spin-0 and -1/2, which simply arises from the ratio of 
squared charges for an isoscalar to an isodoublet state and the u quark dominance. For a spin-1 
target, one would expect the same. However, in the case that (and also TC5) does not contribute 
the factor (3 — 2tG^/Gi)I{2 + (1 — 2TG-i/GiY) gives a dependent suppression. It decreases 
with growing |A^| < 0.5 GeV^, reaches its minimum at |A^| ^ 0.5 GeV^ with ^ —0.1 and then 
increases with growing |A^|. Note the sign change of the beam spin asymmetry, which occurs in 
the region 0.4 GeV^ < -A^ < 0.7 GeVl 

In the following we confront our oversimplified estimates for the ^-integrated beam spin asym- 
metry (jll9|) with preliminary HERMES data from the 2000 run, which have been taken for proton, 
deuteron, and neon target |77] : 

= -0.18 ± 0.03 ± 0.03, (a;B) = 0.12, (-A^) = 0.18 GeV^ (Q^) = 2.5 GeV^ , 
Al^ = -0.22 ± 0.03 ± 0.03, (xr) = 0.09, (-A^) = 0.13 GeV^ (Q^) = 2.2 GeV^ , 
A'l^^ = -0.15 ± 0.03 ± 0.03, (xb) = 0.1, (-A^) = 0.2 GeV^ (Q^) = 2.5 GeV^ , (138) 

where the first (second) error denotes the statistical (systematical) uncertainty. Here the notation 
Al_|_u refers to the fact that this value is obtained from two different data sets of polarized and 
unpolarized deuteron target. Our rather naive estimates give 

Alu = -0.26 , A^^ = -0.34 , A^^ = -0.23 , = -0.21 , (139) 

which are slightly higher than the experimental values. Note that these numbers include uncer- 
tainties induced by the kinematical approximation and, of course, due to the simplification of 
GPDs. Compared to the beam spin asymmetry on a proton target from the 1996/97 run 

Alu = -0.23 ± 0.04 ± 0.03, (xb) = 0.11, (-A^) = 0.27 GeV^ (Q^) = 2.6 GeV^ , (140) 

our naive analytical and numerical estimates, within a certain model that contains also sea 
quarks and twist-three corrections, give in both cases nearly the same value = —0.26 and 
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Alu = —0.27, respectively. The sea quark content gives an enhancement and the exact kinematics 
at twist-tree level a dumping, which results finally in a similar prediction. Nevertheless, these num- 
bers are rather model dependent and in fact by varying the slope -Bsea? ^-parameters, and adding 
dynamical twist-tree contributions by 'hand' our prediction covered the range 0.16 < I^luI ^ 0.37 

For JLAB kinematics with xb = 0.3 [0.2], E = 6 [12] GeV, = -0.25 GeV^ and = 
2.5 GeV^ the quantitative estimates for proton, deuteron, and a isoscalar nuclei A target are 

Alu ~ 0.25 [0.3] , AJ^u ~ 0.2 [0.25], ~ 0.35 [0.45] . (141) 

For the charge asymmetry, there are important differences with respect to the beam spin 
asymmetry. First we recall the non negligible constant term. Moreover, this asymmetry arises 
from the real part of the Compton amplitude and so it is also sensitive to the D-term, which 
drops out in the imaginary part^°. At smaller values of xb the ratio of real to imaginary part is 

TZ = — — — ~ tan (^[2a* — 1 — + pure 'exclusive' contributions/3m (142) 

where the signature Sjr = 1 for the CFFs {H, S, Hi, . . . ,0.5} and Sjr = —1 for the set 
{H, S, Hi, . . . ,Hi}. Here we used that the GPDs behave like H{^) oc S,'"" with a* > 0. Since 
a* ~ 1/2 and a* > 1 for unpolarized valence and sea quarks, respectively, we expect also for larger 
values of xb that the valence and sea quarks give a negative and positive correction, respectively: 

W f < 1 ^ f valence quarks 
R = — — — - = < ^ for <^ ^ . (143) 

ism 7t* I > J 1^ sea quarks 



The D-term contribution has been estimated to be negative, too, and will partially cancel the 
positive sea quark contribution. All these terms enter the Fourier coefficient Aq\ defined in Eq. 
pisp with a reversed sign. The HERMES data 

Ac{(p) = -0.05 ± 0.03(stat) + [0.11 ± 0.04(stat)] cos(0^) , (144) 

for (xb) = 0.12, (-A2) = 0.27 GeV^ (Q^) = 2.7 GeV^ can be explained by two different 
scenarios: 

1. small sea quark contribution and rather small or no D-term 

2. large sea quark contribution and large negative D-term. 

The second scenario predicts a beam spin asymmetry of Alu — —0.37, which is slightly larger as 
the beam spin asymmetry of the 1996/97 run. So a definite conclusion about the D-term can not 
be done with present data. However, we state: if high precision data will be available, the ratio 
of charge asymmetry to beam spin asymmetry allows to pin down the D-term contribution. The 
measured ratio 



4(0) „x 

""O.unp 



0.5 (145) 



^'^A first discussion of tlie imaginary and real part of the Compton amplitude in the context of field theory can 
be found for instance in Ref. 59, 60 . The constant term in the real part, which appears here due to the D-term, 
has been related to a fixed pole in the language of Reggeization. 
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coincides with our estimate p25|) . which yields the value ~ 0.65. 

In both scenarios one expects that the valence quark contribution dominates in the real part 
of the CFF 7i for HERMES kinematics. The real part of this dominant CFF for a proton and a 
scalar nuclei target are 



(146) 



FA(A2) 



where B} and depend on = x-b/{2 — xb), A^, and . As above in the case of the imaginary 
part a similar estimate holds true for a spin-1/2 nucleus target. For larger value of xb, e.g. 
x-Q = 0.12 the asymptotic formula (jl42|) overestimates the size of the real part. In the MRS A' 
parameterization we find i?"™'(^Ar = 0.065) ~ —0.4 and -R°'™'(^Ar = 0.065) ~ —0.25. So we expect, 
e.g., = 2.7 GeV^ the following ratio of charge and single beam asymmetries^^ on a proton 
target 



2~2y + y^ Q^^^^aig^vai 



45'^ i2-y)y 



-0.8 



Q2^«val + QS^'^val 

which is consistent with HERMES data. For an isoscalar nucleus we expect a similar value: 



(147) 



A, 



A(l) 



2 - 2?/ + /?«valg«val + /^rfvalgrfval ^ . Q2^i 



A 



A(l) 
LU 



y)y 



0.7. 



148) 



However, these ratios should be smaller in reality, since the sea quarks enter the numerator and 
denominator with negative and positive sign, respectively. 

The ratio of charge to beam spin asymmetry for a spin-1 target has a rather complex model 
dependence. Even if we set 7i™^ and to zero, we should worry about the unknown 'quadrupole' 
contribution TC3, induced by the sea quarks: 



A 



A(l) 



A- 



A(l) 
LU 



^ ^ 9 R«val^«val _|_ Prfval^rfval _|_ O '-'I 

2-2y + y^^ 1 ^ +^Gi 


\ 2TGf='(Gi-2rG3)" 
G|<=='{3Gi-2rG3) 


Dsea 2 JZi Qi 9' 


(2 - y)y QUval + O'^val + 2^ 

^ ^ Gi 


" 2TGf-(Gi-2rG3)" 
G|<=='{3Gi-2tG3) 


/-^ i — u,d ^ i 



(149) 



Here we have set for simplicity T-Cf^/T-Cf^ = Gf^/Gf^ = 1 — + Q^£^ and assumed the same 
^sea g^jj ggg^ quark species. If the sea quark contributions are small, we would have the same 
ratio ()148|1 as for isoscalar nucleus with spin-0 or -1/2. 

It is interesting that for larger values of xb the sign of the real part of the CFFs might be 
changed. In our model that happens for valence u[d\ quarks at xb ~ 0.3[0.2], while the sign 
of sea quarks and D-term remains the same. There is also a significant difference between these 
contributions. While the real part of the sea quarks is expected to drop rather fast with increasing 

^^This ratio has the advantage that the BH contribution and partly the kinematical prefactors of the interference 
term drop out and so it is after restoration of the y dependence exact up to the neglected squared DVCS term and 
the approximation of the CFFs. H and E CFFs, which are kinematically suppressed, are still neglected. 
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xb, the D-term gives a constant contribution. JLAB kinematics would be suitable to explore this 
region and provide so more insight in the structure of GPDs, unfortunately, no positron beam is 
available. 

Now we come to the longitudinally polarized target spin asymmetry. For a proton target the 
ratio of this target spin to the beam spin asymmetry is 



^ul(0) 2-2y + y' 
Alu(0) ^ i2-y)y 



2 Fi 



(150) 



Both terms in the squared bracket on the r.h.s. are positive and become smaller^^ with decreasing 
xb- So within our simplifications we expect the same sign as for the beam spin asymmetry and a 
rather small value for not so large xb or small y. In the case of a isoscalar nuclei spin-1/2 target 



2-2y + y 



AU<P) {2-y)y 



g«val -|- q 



+ 



xbFI 
2AF( 
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the form factor dependent part will be suppressed by the atomic mass number A. As long as this 
term is not enhanced due to a node in F^(Aq) at given Aq, the ratio of polarized to unpolarized 
quark distribution dominates and we expect for xb ~ 0.2, Q'^ ^ 2 GeV^ and HERMES [JLAB@12 
GeV] kinematics 



2 - 2y + ?/2 Ag"vai + Ag'^v, 



0.5 [0.2] . 



At^i^) {2-y)y g«va. + g^va. 

If we neglect 7^3 and Ti^ CFFs, we find for an isoscalar spin-1 target a similar prediction: 
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which differs from the previous cases by a A^ dependent prefactor. Unfortunately, this ratio is 
also rather sensitive to the details of and GPDs, compare Eqs. ()121|) and ()127|) . and so no 
definite estimate can be given. 

Finally, we discuss the tensor polarization asymmetries. In the case that Ti^ and 7^5, are small 
we find that these asymmetries can be expressed by the charge and beam spin asymmetries: 



AtA^) A 



A(l) 



2 {-T + AiT-?>e)G^ 

3 G\ — 7;TG^ 



(154) 



For small values of — A^ and ^ these ratios are rather small, e.g., — A^ = 0.2[0.1] GeV^ and 
Xb = 0.1 we find the number ~ 0.2[0.08]. A significant deviation from this prediction would 
indicate that Ti^ and/or Tis are comparable in size to Tii. For instance, if would be small and 
assuming that Ti^ = G^/GiHi this ratio is roughly double so large. 
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The ratio of form factors is for small value of |A^| given by the anomalous magnetic moment F2/F1 w 1.793(1 — 
A^/1.26), while the ratio of polarized to unpolarized quarks is estimated from the measured ratio of polarized 



to unpolarized structure functions in deep inelastic scattering which behave like 
gi/Fi ~ 0.2 for xb ~ 0.1 . 



XB for increasing xb, e.g. 
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all targets 


proton 


neon 


deuteron 


model 


Kal 


^sca 


5sea [GeV-^] 


'^sca 


Ssea [GeV-^] 


fisea [GeV-^] 


A 


1 


OO 


9 





9 


20 


B 


OO 


OO 


9 


-3 


12 


30 


C 


1 


1 


5 





6 


15 



Table 2: Parameter sets for H and Hi GPDs of a spin-0 nucleus, proton, and deuteron. 
5.2 Numerical results 

In the previous section we tried to convince the reader that a rather simple GPD ansatz explains the 
measured DVCS asymmetries. To avoid a misunderstanding, we do not claim that our estimates 
are precise. Rather we expect that more realistic GPD models yield large deviations from these 
naive estimates. Nevertheless we like to give an impression of this model dependence in the 
following. We also have to keep in mind that the kinematical approximations are rather rough. In 
the following we correct this and take the exact expressions for the BH cross section and K factor, 
defined in Eq. fl23j) and the full ^ and dependences in the twist-two predictions of interference 
and squared DVCS term. 

First we consider the beam spin asymmetry measured at HERMES. For the proton target 
we rely on the GPD models A, B, and C, employed in Ref. The parameters of these three 
models are listed in Tab. |21 Compared to model A, the sea quarks are suppressed and enhanced 
in model B and C, respectively, by a larger slope -Bsea and 6sea oo. We also suppress in model 
B the valence quarks by taking 6vai oo. At twist-three level the model A and C are based on 
the WW approximation, while in the model B quark-gluon-quark GPDs are added by 'hand', see 
Ref. for details. We neglect the D-teim, which affects the observable in question only sligthly 
at twist-three level. Taking the mean values of the 2000 data f|138|) . the ^-integrated beam spin 
asymmetry ()119j) at twist-two and twist-three level reads: 



[-0.29,-0.25,-0.41 ] ^ [twist-two] ^ rAT.r.i 
^Lu = S W for < , > and models A, B, C . 

\ -0.29,-0.20,-0.37 J \ twist-three J L > > J 



(155) 



As we realize the model B predictions are consistent with the experimental data, while the model 
A gives a slightly too high value for I^luI- We remind that both model predictions are consistent 
with the 96/97 data. The large value of the model C prediction is caused by the sea quarks. 

Now we present the predictions for the neon target at twist-two and -three level. Here again 
the sea quark contributions are relatively suppressed in model B and C, see Tab. |21 and we again 
add for the B model a quark-gluon-quark GPD contribution by 'hand', in the same way as it has 
been done for the proton. We find with these models: 



(156) 



f 1-0.30. -0.28. -0.54n f tw,st-two 1 |a, b, C] . 

\ [-0.34, -0.24, -0.49] J \ twist-three J L ' > J 

Again the model B prediction is on the one a level consistent with the experimental data 

Now we discuss estimates for an unpolarized and longitudinally polarized deuterium target in 
more detail. The models for the reduced deuteron GPD Hi are listed in Tab.|21and they share the 
same qualitative features as described above. The D-term contributions are neglected and we set 
for simphcity /igea = 1 and Qsca = 0. To illustrate possible bound state effects, we sometimes take 
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Figure 4: Model dependent estimates for the cross section, integrated over the azimuthal angle 
(j), for the scattering of a positron [electron] beam on a deuteron target at Eg = 27.6 [12] GeV, 
Q2 = 2.5 [2] GeV^ = -0.1 GeV^ (a) [(c)] and = -0.3 GeV^ (b) [(d)]. The squared BH 
contribution is displayed as bold dotted line, the DVCS cross section as dash-dotted (dotted) and 
the leptoproduction cross section for positrons as solid (dashed) curves for model B (model C). 

into account, equating the reduced GPD with that of Hi, or alternatively H^, equating it with 
the antisymmetric part of Hi{x, rj, A^). These two additional sets will be denoted with a prime or 
a hat, respectively. As mentioned above, we neglect the kinematically suppressed contributions 
H2, H4, H2, and H3. The model of Hi is based on the GS A parameterization |^ with bsea, = 1. 
If not stated otherwise, we equate the slope parameters 5vai = -Bsea = 20 GeV~^, which is roughly 
the size of the slope parameter for the charge form factor Gq. 

Let us first consider the size of the unpolarized cross sections for the positron and electron 
scattering off a deuteron target at HERMES and JLAB kinematics, respectively. It is demon- 
strated in Fig. Elthat at larger value of xb, i.e., smaller value of y, the cross section is dominated 
by the DVCS one, however, its size is comparably small. With decreasing xb, the BH cross section 
starts to dominate over the model dependent DVCS one. In model C the sea quarks also induce 
a growing DVCS cross section due to an associated large imaginary part in the DVCS amplitude. 
Typically, the ratio of BH to DVCS cross section and the total cross section is larger at JLAB 
compared to HERMES kinematics. 

In Fig. El we show the model dependence for the predictions of the beam spin asymmetry 
Ai^±. In panel (a) we confront them with the preliminary HERMES measurement ^T7\. However, 
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Figure 5: The beam spin asymmetry Al±(0) (a), Al± as function of xb (b), — (c), and (d) 
[without evolution] for the scattering of a positron on a deuteron target at Ee = 27.6 GeV. The 
fixed kinematical variables are xb = 0.1, = 2.5 GeV^, and = —0.2 GeV^, where the models 
A (dash-dotted), B (dashed), B (solid), and B' (dotted) are specified in the text. 



we remind that the analysis is based on two different data sets that, certainly, also contains 
incoherent scattering events. Model A (dash-dotted) gives a larger asymmetry than model B 
(dashed), caused by the larger sea quark content of the former model. Note that for the same 
reason, model C predicts an even larger asymmetry, which is not displayed. If we take into account 
the H3 contribution via the B' model (dotted), the beam spin asymmetry in panels (a,b,d) is of 
the same size as the model A predictions. It is displayed in panel (c) that model A and B' are 
distinguishable due to the A^ dependence. In absence of i^s we observe for the typical HERMES 
kinematics only a rather small difference between the asymmetries and Ai^±. The latter 
asymmetry is sensitive to and we find a large reduction for the model B (solid). We recall that 
the 0-integrated asymmetries, displayed in panels (b)-(d), are smaller than the non-integrated 
asymmetries at = 7r/2. This integration also changes the A^ dependence^^. For JLAB kinematics 
we find the same qualitative features, where of course for electrons the sign of the asymmetry is 
reversed. Let us add that within the model B we have a slightly larger value than in our result 

^^The unintegrated asymmetries H112(l and H115|l are roughly spoken proportional to y^—A^/ and thus will 
constantly decrease with decreasing |A^|. For the integrated asymmetry we rather observe a growing that comes 
from the (j) integration, indicated in Eq. (|119|1 . Of course, in any case the asymmetry vanishes in the limit 
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Figure 6: The longitudinal target asymmetry Aui^{(f)) (a,c) and Aul (b,d) versus (p and xb, 
respectively. For HERMES (a,b) [JLAB@12 GeV (c,d)] kinematics we take xb = 0.1 [xb = 0.2], 
-A2 = 0.2 GeV^and = 2.5 GeV^ as well as employ the model A with 5^ai = |sca = 20 GeV"^ 
(dash-dotted), E^ai = 30 GeVj^ B,^^ = 10 GeV'^ (dashed), S^ai = 10 GeV~^ B,^^ = 30 GeV'^ 
(dotted), and B with E^ai = 5sca = 20 GeV"^ (solid). 



(jl39|) of the previous Section, caused by both kinematical and sea quark contributions. The other 
models induce the typical effects, we have already discussed: 

(y)( 1-0^37. -0,26. -0,23, -0.34] 1 . 

l^LiJ [[-0.34,-0.24,-0.16,-0.34] J L . . . J v ; 

These results will also be altered by twist-three contributions, i.e., change of the normalization 
and excitation of higher harmonics. A quantitative estimate of such contributions is beyond the 
scope of this paper. One might expect a similar change as in the case of a proton target, discussed 
in great detail in Ref. [7j. 

As discussed above, see also Eq. (I153|) . the size of the longitudinal target spin asymmetry Aul is 
rather model dependent and allows to access the imaginary part of the parity odd OFF Tii. This is 
demonstrated for HERMES and JLAB kinematics in the panels (a,b) and (c,d), respectively, of Fig. 
ini Here we set xb = 0.1 and xb = 0.2, respectively, as well^as — = 0.2 GeV^ and = 2.5 GeV^. 
For models A (dash-dotted) and B (solid) with same Hi, taking Svai = -Bsea = 20 GeV~^, we 
find a small asymmetry. Within this model parameter the contribution of Tii is negligible, since 
its valence and sea quark contributions partly cancel each other. Thus, we observe in agreement 
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Figure 7: Tensor asymmetries Azz{4>) (a,c) and Ai^zz{<P) (b,d) for HERMES (a,b) and JLAB@12 
GeV (c,d) kinematics within the same GPD models as in Fig. Eland mean values as in Fig. IHl The 
thick dotted line in panels (a,c) shows only the BH cross section. 



with Eq. ()153p that the asymmetry has the same sign as the beam spin asymmetry and generally 
increases with growing xb, as it is displayed in panels (b) and (c). This is mainly caused by the 
fact that the asymmetry is proportional to 1/y. Thus, comparing panel (a) and (c), we realize 
that for the typical JLAB kinematics (larger xb) our estimates for |yluL(0 = 7r/2)| are slightly 
larger than for the HERMES one. If the polarized valence quarks are relatively suppressed with 
respect to the sea quarks, displayed by the dashed line, the negative sea quark contribution takes 
over and induces the sign change in panel (a). For larger value of x-q, see panel (b), the asymmetry 
turns out to be positive. For the contrary case that the valence quark contribution is enhanced 
(dotted line), the target spin asymmetry becomes sizably large. 

The sign of the charge asymmetry Ac is not predictable. For the mean values used above, 
we find that the model A gives an asymmetry that is negligibly small, the B models provide a 
positive and the C model a large negative contribution, where no D-term has been taken into 
account. This observation is in qualitative agreement with our results for the proton target and 
affirms once more that an access to the D-term contribution is intricate. 

The tensor asymmetry Azz{(f)), cf. Eq. ()116|) . of an unpolarized positron and electron beam is 
displayed in Fig. [7|^a) and (c) for the same HERMES and JLAB@12 GeV kinematics as in Fig. IHl 
respectively. The thick dotted line shows the pure BH contribution, which is in correspondence 
with Eq. (jl28j) almost fiat for HERMES kinematics (a). In the case of JLAB kinematics (b) its 
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contribution is much smaller compared to the approximation fll28p . which is caused by higher 
order terms in ^ and r that have not been taken into account [see analogous discussion of Eq. 
P53p ]. Subtraction of this contribution gives a dominant cos(0) contribution that is proportional 
to the real part of the linear combination (jl29|) . As in the case of the charge asymmetry the sign 
depends on the details of GPDs. In panels (b,d) we have ploted the tensor asymmetry Al2z(0), 
cf. Eq. ()131|) . for a polarized lepton beam, which is essentially given by a sin(0) harmonic. Its 
amplitude is proportional to the imaginary part of CFFs combination that enters Eq. ()132|) . If 
both CFFs 7^3 and 7^5 are small, as it is realized in models A (dash-dotted) and B (dashed), 
we find only a small asymmetry, which has the same sign as the beam spin asymmetries. This 
is in line with approximation ()154|1 . This tensor asymmetry is rather sensitive to bound state 
effects. For instance, in the case that there is an important QmT-C^ contribution the asymmetry 
becomes sizable. On the other hand, if we assume a large positive QmTi.^ and a negligible QmJ-Cs 
contribution, i.e., taking model B' (solid), we find a sizable asymmetry that has a reversed sign. 

6 Summary 

In this paper we considered the leptoproduction of a photon on nuclei up to spin-1. For spin-0 
and spin-1/2 targets the theoretical predictions in terms of nuclei GPDs simply follow from the 
known results, presented in Refs. [Ul ESI U\ at the twist-three level, by appropriate replacements 
of form factors, GPDs, and kinematical variables. As a new result we added to this collection the 
azimuthal angular dependence of the leptoproduction cross section of a photon on a spin-1 target 
at twist-two level to leading order in perturbation theory. Since the hard-scattering part is for all 
targets unique, the result for the considered harmonics can immediately be extended to NLO. At 
this order also a gluon tensor contribution appears, which induces a cos (30) and cos (20) harmonic 
in the interference and squared DVCS term, respectively. We should add that the elaboration of 
the twist-three sector, providing the zero and second harmonics in the interference term and the 
first ones in the DVCS cross section, is straightforward in our covariant formalism. 

For a spin-1 target nine CFFs enter the twist-two sector, which are given as a convolution with 
nine GPDs. From the theoretical point of view, the measurement of the imaginary and real part 
of all these CFFs is possible for the case of a polarized lepton beam and target, with an adjustable 
quantization direction, for both kinds of leptons is available. Moreover, an appropriate Fourier 
analysis allows to eliminate the twist-three contamination. In this way one has the maximal access 
to the twist-two deuteron GPDs, given at the diagonal x = ^ and as convolution. 

An important issue is the contamination of the leading twist-two prediction by power sup- 
pressed contributions. Naively, one would expect that such contributions scale with A^/ and 
Q^. Fortunately, within our model we showed that the so-called target mass corrections have 
the same A scaling as the leading twist-two contributions. A more detailed investigation of power 
suppressed contributions would be very desired. 

Certainly, a dedicated experiment is required to measure all CFFs. Having on hand only an 
unpolarized target one can access the beam spin and charge asymmetry via the Fourier coefficients 
■^f.unp '^unp' given as a certain linear combinations of Hi, Hs, and TC^ CFFs. A longitudinally 
polarized target allows to measure four further independent observables. Two of them, e.g., the 
tensor polarization combined with beam spin or charge asymmetry, are predicted by another linear 
combination of the afore mentioned three CFFs, entering the Fourier coefficients and llp- 
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Unfortunately, such experiments do not aUow to access all three GPDs in question. The remaining 
two observables are related to the target spin flip odd asymmetries, given by the FCs Lp and 
Lp. In a single target spin flip experiment one can measure Lp and get access to the imaginary 

part of Til- On the other hand, combining the data for polarized target with charge asymmetry 
or the flip of the beam helicity allow then to measure Lp and so also the real part of Tii. In all 
these observables, we discussed so far, the remaining flve CFFs are kinematically suppressed or 
do not enter the theoretical predictions. 

To estimate the size of different observables, we proposed a simple model for nuclei GPDs that 
is based on the A-scaling of the longitudinal degrees of freedom and a factorized ansatz of 
the nucleon GPDs. For a spin-0 or -1/2 nucleus target, this ansatz is sufficient to fix all reduced 
GPDs in terms of the proton ones (of course, they are also unknown). The remaining degree of 
freedom concerns mainly the form factors, especially, in the sea quark sector. Of course, binding 
effects should alter this oversimplified ansatz, e.g., we expect that t-channel exchange forces can 
induce an important contribution in the 'exclusive' region of GPDs. Such effects might be visible 
by a precise measurement of single spin and charge asymmetries. For a spin-1 target the large 
number of GPDs opens a new window to study nuclear binding effects. Compared to deep inelastic 
scattering or elastic lepton-deuteron scattering, they provide additional information, contained in 
the CFFs. Beside the modifaction of the the scaling law for the CFFs Tii and Tii, they induce 
no n- vanishing CFFs H3 and H5. We conclude that nucleus GPDs give a new testing ground 
for non-perturbative methods that have been used in nuclear physics for the calculation of form 
factors, see reviews [301 Ell- 

Finally, we estimated for the typical kinematics of present fixed target experiments, different 
asymmetries and showed that a qualitative understanding of them is possible by means of ana- 
lytical formulae, which are obtained by a simplification of kinematics and GPD models. For a 
deuteron target such an approximation allows to discuss the contribution of the CFFs Tis and 
7-^5. The LO analyses of the pioneering measurements of DVCS [12112111221 on the proton suggest 
for — ~ 0.3 GeV^ the dominance of valence-quark GPDs in the valence quark region with no 
essential enhancement by the skewedness effect. Assuming the same features also for the nuclei 
GPDs, our estimates agree with the ^-integrated beam spin asymmetries measured on neon and 
deuteron targets. In the latter case this asymmetry is already sensitive to the CFFs TC3 and H5. A 
measurement of the tensor polarization would certainly provide a new constraint for these CFFs. 
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Sections El and El 



A Parameterization of the electromagnetic deuteron form 
factors 

The set of form factors Gi, G2, and may be expressed by the charge monopole, magnetic 
dipole, and charge quadrupole form factors: 
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Gq(A2) = Gi(A2) - G^iA') + (1 - r)G3(A2), 
The normahzations of these form factors read 
Gc(0) = 1, GMiO) = /i, = 1.714, Gq(0) = = 25.83. 
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Their parameterization, which is inspired by the counting rules for large — A^ was taken from 
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Here we use the standard dipole parameterization for the nucleon form factor 
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0.71 GeV^ 

and the helicity transition amplitudes in the infinite momentum frame are given by 
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Counting rules predict for large — A^ the following behavior 
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Table 3: Fitting parameter sets for the electromagnetic Deuteron form factors. 

To reduce this set to twelve parameters, one may introduce for each group a,, and 7i the 
algebraic relations: 

i — 1 

a\ = \ {a\ - a^) — — for i = 1, . . . , 4. (164) 
The fitting parameters are taken from Ref. [S3] and are given in Tab. El 

B Results for the Fourier coefficients 

Below the twist-two results for unpolarized and longitudinally polarized target are listed. All 
results have been expanded for small r = ^j^, i.e., for <C M^. Terms proportional to 7^3, G3 
are given up to order 0{t\ because is at = roughly 20 times larger than the other 
form factors. The BH amplitude squared has been exactly calculated. However, the results are 
very cumbersome and, thus, they are only presented in leading order of 1/Q^. We note that the 
rather lengthy results for A^o'lu^^ ^'^'^ ^(?llp '^^'^ ^e obtained through the much shorter results 
for unpolarized and transverse-transverse coefficients via equations (jUj). 

B.l Interference term 

Here we give the matrices A^f ^np; -^f,LP; ^^'^ ''^f.TTPE; which appear in Eq. 
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B.2 DVCS amplitude squared 

The matrices -A/^°7^^ entering Eq. (gHD, for i G {unp, TP-, LLP, LTP+, TTPS, TTPA} are of 
the form 







rDVCS 



with Ajj^y*-'^ being a 5 x 5 and BqJ'-^^ being a 4 x 4 matrix. The matrices A^Jj^y^^^ for j e 
{LP, TP+, TTP±} have the building blocks 
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where C^J*^^ is a 5 x 4 and D^J"^^ is a 4 x 5 matrix. 

We list only the sub matrices A^J^^, Bjj'J^^, C^J^^, DqJ'^^ that are needed for an unpolarized 
and longitudinally polarized target: 
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C Estimate of target mass corrections 



Target mass corrections for a scalar and spin-1/2 target have been evaluated in Ref. [SZI, where 
A^/ corrections have been neglected. This result will serve us to derive the scaling law of the 
target mass corrections with respect to the atomic mass number A for a general Compton process 
in the light-cone dominated region. It is sufficient to consider the simplest case of a scalar target. 
Then the target mass corrections in the parity even sector read for the Compton amplitude 



^liiiV) = I dydzh{y,z) 



+0 



J [z^ -z 



An analogous formula for the other leading twist amplitude J^2 holds true. Here h{y, z) denotes 
the DD and the coefficients 

Cf\s-') = -Y^. and C«(--) = ^^^ + 21n(^) , (166) 

depend on the variable S = ^/{y + rjz). Let us set the mass M and the scaling variables in Eq. 
(jl65|) equal to the nucleon ones: Mjsij^jsiyVN- The Compton amplitude for a nucleus A at smaller 
values of the Bjorken variable^^ follows from the replacements: 

Mn ^ M = AMn , ^ ^ ^n/A, r]N V ^ Vn/A, and h ^ . (167) 

Corresponding to our nuclei GPD model (f7n|) and the GPD definition ()50|) in terms of a DD, we 
read off the scaling law for the nucleus DD: 

h'^iy, z) oc e{\Ay ±z\< \)h{Ay, z) . (168) 

Changing the integration variable y = y'/A and renaming y' — > y, we find that the Compton 
amplitude for a nucleus target reads 



Ff{^,v) oc / dydzh{y,z) 



\ J Q\y + VNzY V iN 



J [z^-z 



(169) 



As we realize by comparison with Eq. (|165p . setting M = M^, ^ = ^n, and r] = ri]\f, the target 
mass corrections do not scale with A in respect to the leading twist-two term. It can be shown in 
the same way from Eqs. (30) and (31) in Ref. [HZ] that also the resummed target mass corrections 
possess the same property. From the same representation it also follows that this statement is 
true for a target with non-zero spin content, including the parity odd sector. 

Beside this kind of corrections also dynamical ones, which are expressed in terms of multi- 
parton correlation functions, appear in the CFFs. Unfortunately, only little is known about them 
and it remains an open task to consider the full twist-four sector. 

^"^We only use this condition to simplify the discussion. For larger values of xb one should take the relation 

e = (i+e)ejv/(i + eiv)A 
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